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Eine früher abgeleitete algebraische Gleichung ($ 1) für den gegabelten Verdichtungsstoß wird auf Stöße 
in Luft angewandt. Neben numerischen Daten für besonders interessierende Stoßkonfigurationen, die als 
spezielle Punkte und Kurven im W eise-Diagramm für Luft auftreten, werden deren Gleichungen für beliebige 
x aufgestellt. ($ 2—4). In $ 5 findet man schließlich eine kurze algebraische Behandlung der Herzkurven- 
methode und der Schnittzahl von Haupt- und Nebenherz. 

An earlier derived algebraic equation for the bifurcated compression shock ($ 1) is applied to shocks in 
air. Besides numerical data for especially interesting shock configurations, which are represented by special 
points and curves in he Weise-diagram for air , their equations for arbitraryx» are derived ($ 2—4). 
Finally $ 5 contains a short algebraictreaiment of ihe heart curve method and of the number of intersec- 
tions of a principal and a secondary heart curve. 

Une Equation algebrique, derivee plus 1öt pour les chocs de compression bifurgues ($ 1), est appliquee 
aux chocs en air. Aux configurations particulierement interessantes, qui sont representees par des points ou 
courbes sp&ciaux du diagramme de Weise, des valeurs nume£riques sont donnees pour x de l’a ir ; mais leurs 
Equations sont derivees pour chaque x($2—4). Enfin on trouve en $5 un breftraitementulge- 
brique de la methode des courbes en coeur et du nombre des points d’intersection des coeurs principaux 
ei secondaires. 


BeiBefeHHoe npeikze anreöpanmyeckoe ypasHeHnue ($ 1) Ma YINOTHAMOINero ÖndypKaumoH- 
HOTO yAapa IIPHMeHAeTca K YAapaM B BOsAyXxe. Hapaıy c UNCJIOBBIMH MAHHBIMH IH OCOÖeHHO 
HHTePecHbIX YAAPHBIX KOHPUTYaIMa, BOSHUKAMNIHX B BUNE OCOÖBIX TOYeK HU KPHBBIX B HA- 
rpamMe Beüse A1A BO3AYyX&a, YCTAHABAIHBAWTCH UX ypasHeHHun Ana NWOBIX x ($2—4), 85 
cOAEeP:KUT, HAKOHEN, EPAaTKyO ANTeÖPAaHuH4UecKyI TPAKTOBKY MeTola KapAHoufsI MH yucıa 
ceyeHMÄ TIABHOÄH HU BTOPOCTEIIEHHOH KAPAHOHN. 


$ 1. Eine algebraische Gleichung für den Gabelstoß 


In einer früheren Abhandlung [1] wurde aus den 2 Gabelbedingungen von A. Weise[2] 
6=6, +9; und 7Y=IL,: I eine algebraische Gleichung 6. Grades für die gegabelten Ver- 
dichtungsstöße in stationären Überschallströmungen abgeleitet. Die Koeffizienten A, dieser 
Gleichung hängen von 3 be- 
kannten Variablen ab. Die Hauptstoss 
erste, k—= (x — 1)/(x +1) mit (h 
* = Cp/6,, beschreibt die Gas- 
art. Die zweite, 1=ctg? a—= 


M2 —1 = (# —1)/(1—ko}), ist 
ein Maß für die dimensions- KA, 153 
lose Anströmgeschwindigkeit Stromlinie (0) 


@, —= w,/a* vor dem Haupt- 
und dem vorderen Nebenstoß. 


P, __-<Trennungs- 


&ist der Machsche Winkel, nströmung ——> linie 
M die Machsche Zahl in der Ans g SER 
Anströmung. Die dritte, m= w,p, 9,7, S 


hängt von der dimensions- 
losen Normalengeschwindigkeit 
On, = Wn,/anin gleicher Weise 
ab wie 1 von o,. Sie ist ein Maß 
für die Neigung des Haupt- 
stoßes zur Anströmrichtung. 
Dies zeigt ein Blick auf die For- 


: 
mel tg? 9,—(1+ m)/(I— m) für Pza, 92a,'za 


den Winkel #, zwischen den bei- (a) Nebenstösse I(b) 

den. Die Unbekannte der Glei- vorderer hinterer 

chung schließlich = M,,— 1 

Tun = n i i ichtungsstoß. eschwindigkeiten, 
EEE EB Begsienungen su: Gagsbelten Tertentunmetos;  üessheintieketten, 


zu l und m analog aus der mung und Stoß, IT Quotient des Druckes nach dem Stoß durch den Druck vordemStoß. 
6 
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dimensionslosen Normalengeschwindigkeit @,, , = Wnı ,/@r, „ für den vorderen Nebenstoß ge- 


bildet. In tg? 9,. = (1 + @)/(l— x) beschreibt sie die Neigung des vorderen Nebenstoßes zur 
Anströmrichtung. a*ist die kritische Schallgeschwindigkeit, a, und a,;, sind die kritischen 


Normalengeschwindigkeiten für den Hauptstoß und den vorderen Nebenstoß. Weitere Be- 
zeichnungen gehen aus Bild 1 hervor und eine ausführliche Formelsammlung findet man in 
I1$2 und in IL[3]. 
Wir definieren drei Polynome 
A2)= — ka — a, + a,2+ (I— m), Ba)=br+bs2+b, Cl)= 5b — 18 — co (1) 
mit den Koeffizienten 
a4 —=—14+(1+k)l—2kmg =1+k—kl+ km,b,—= (1+ m) (I— m) + (1+ km), 
b,=(1+k)b,— mb, db, = k(1+1) — (1 — k?)m, = (1+ 1) (I — m) m, 
G=(—m)[1+1+(—1+2%k+lm]. 
Damit läßt sich obenerwähnte Gl. 6. Grades in der Form schreiben: 
{(b,©+ b1) C(@)+ bo(bo® — 1) — (1 +1) b[A(a)+ m(k?+ 2 — a)]}? (2) 
— 4(14+ m) [1 +1— (1 — k)m]? b,(1— x) A(x) B(x) en 
Diese Gl. ist dieselbe wie Gl. (113) in berichtigter Form [4]. Ordnet man sie nach Potenzen von 
x, so erhält man 
(8,2348, +82+8) = R(T® + T,#? + T,#+ T,0+ T,&®+T,c+T,) : - - 8) 
mit den zwölf von k, I, m abhängigen Koeffizienten 


S—=b[l&»+k(l+D)), = buldi+ (+1) (aa — km)] — bzcı, 

Ss = bl — (+) (+ wm)] — ba — bc, SH = — bla + (1+ DM (l— m — ma,)]) — co, 
R=4b, 1+m)[1+1— (1— km], T,= kb, T,=kb, +b,(,— kl), 

T,— kbo+ b(a, — kl) — bz(a, + 1a,), Tz— bu(a — kl) — b, (a, +la,) —b,(ı— m — la,), 
T,— — bu(a, + la,) —b,(I— m —la)+lli— mb, T,=—bu(l — m — la,)+ I(l— m)b, 
T,—= (1 — m) by. 


Vereinigt man die Glieder gleichen Grads in & auf beiden Seiten von (3), so entsteht die Gl. 


6 
DER I a Bi NE 
n=0 
mit 


GEBIRT,G=2 88, RT, Sr 28,8, RT OS SS ran 
.2=28%8+ 8 —RT, 0, =28 5 —RT,, = S8—RT,. 


Diese Koeffizienten (', sind gemeinsame Vielfache der Koeffizienten A, von Gl. (114) 


6 
> A, FE — 0 ©. 0, er Did. zul Ve we u ae Ze (5), 
n=0 
nämlich 0, —= m? A,„,n=0,1,...6. Eine explizite Darstellung der Koeffizienten A, findet man 
in I hinter Gl. (I14). Physikalisch brauchbar sind solche Lösungen x von (5), für die 
0<zr<m<siögilt. 


Anwendung der Gleichung auf Luft (x = 1,405). 


Man unterscheidet zwei Lösungsbereiche, die wir in Bild 2 als I und II darstellen. Wir 
vereinigen beide Bereiche dadurch zueinem Bild, daß wir Abszisse und Ordinate für Be- 
reich II gegenüber Bereich I vertauschen. Für Bereich list x auf der Abszisse, m auf der Ordi- 
nate abgetragen, für Bereich II umgekehrt. Beide Lösungsbereiche hängen an der 45°-Linie 
© — m zusammen, welche die Bedingung für den Grenzfall ‚‚hinterer Nebenstoß ist Mac h sche 
Linie“ bedeutet. Sie überschreiten in dieser Darstellung nicht die 45°-Linie, eine mehrfache 
Überdeckung findet nicht statt. Die ausgezogenen Kurven sind Kurven konstanter Anström- 
geschwindigkeit ®, bzw. 1=const. Man entnimmt aus ihnen die m- und z-Werte für einen 
Kurvenpunkt, trägt sie in die Formeln des gegabelten Verdichtungsstoßes ein und erhält ein 
physikalisch mögliches Strömungsbild. Die gestrichelten Kurven in Bild 2 haben folgende Be- 
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deutung: Punkte auf (, führen zu einer Gabel mit senkrechtem Hauptstoß, Punkte auf (0, zu 
einer mit senkrechtem hinterem Nebenstoß, für Punkte auf C, wird im Hauptstoß maximale 
Ablenkung von gleichem Betrag wie im 
einfachen schiefen Stoß (für das be- 
trachtete ®,) erreicht, die Kurve (©, 
bzw. (©, verbindet die Maxima in x- 
bzw. in m-Richtung des 2. Lösungs- 
bereichs, über sie gibt $5c bzw. $2f 
weiteren Aufschluß. Die Bedeutung der 
besonders markierten Punkte ist in allen 
Bildern gleich (dieselbe Stoßkonfigura- 
tion!). 


$ 2. Der zweite Lösungsbereich für 
o,< 1,85 

a) Im Gebiet der gleichdrehenden 
Vorwärtsgabeln (,>0) für Luft stim- 
men die Ergebnisse der verschiedenen 
Autoren über das Aussehen des ö, ö,-Dia- 
gramms (Bild 3) im Bereich der Anström- 
geschwindigkeiten 1,672 < o, < 1,84172 
nicht überein. A. Weise [2] hat 
diesen Teilbereich offengelassen !). W. 
Wuest [5] war der Meinung, daß die 
Kurve für &, = 1,84172 ohne Unter- 
brechung in den Ursprung einlaufe und 
daß es einen (in seinem Bild 3 schraf- 
fierten) Bereich gäbe, in dem keine 
Lösung existiere. Fr. Wecken [6] 


K 
@,=1,84172 


R . Bild 2. Kurven des 1. 1 2. Lös sb ichs (I und II) für Luft bei 
hat diese Aussage Wuests wider- ild urven des 1. und wen ereichs (I un ) für Lu ei 
18° T or 


-&-Gerade\ 


7 Der Fer SE 


o£ iS 
f 
Bild 3. 6, ö,3- Diagramm (Weise-Diagramm) des 2. Lösungsbereichs für ®, < 1,85. Die Zahlenangaben an einzelnen Kurvenpunkten 
bedeuten die zugehörigen m-Werte in Bild 2. Die gestrichelten Kurven verbinden gleiche m-Werte auf verschiedenen 
BR N IIRR-: ©, = const-Kurven 
1) A. Weise, Jahrb. d. deutsch. Luftfahrtforsch. 1943 IIA 011, S. 8, Bild 10. An 2 Kurven des er- 
wähnten Bereichs (rechts der ö = ö„-Geraden) muß die Beschriftung geändert werden: statt 1,8 bzw. 2,0 heißt 
es 2,0 bzw. 2,2. Dies geht auch aus dem zugehörigen Text auf S. 8 links hervor, wenn man dort in der 2. Zeile 
von oben einen Druckfehler in », = 1,84172 (statt 1,184172) beseitigt. Der richtige Wert steht in der letzten 
Spalte von Tafel2. Zur Frage der Zweischnittigkeit vgl. jedoch Fr. W ec ken[6]oder $ 5b dieser Abhandlung. 


6* 
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legt. Zur Nachprüfung der Sachlage und als Anwendung unserer (Gl. 4) oder (5) haben 
wir diesen Teilbereich nochmals untersucht. Die Rechnungen wurden mit x = 1,405 
d.h. mit k = 0,168 399 17 unter Benützung einer Rechenmaschine Mercedes Euklid 38 M 
durchgeführt. Wir bestätigen in den wesentlichen Punkten die Aussagen Fr. Weckens. 
Bild 3 gibt unsere Ergebnisse wieder. Um einen’ Doppelpunkt E für (o,)z = 1,672 333 2 auf 
der ö-Achse des 6, ö,-Diagramms bei öz — 12° 46’ 51” (statt 19,4’ nach Wecken[6], Tafel 2) 
herum gruppiert sich eine Folge von ®, = const-Kurven mit wachsendem @,, deren „Fuß- 
punkte“ auf der ö-Achse stetig auseinanderrücken und deren Maxima immer höhere ö,-Werte 
erreichen. Über die Daten des Punktes ZE und der Fußpunkte der Nachbarkurven mit etwas 
höherem , gibt Tabelle 1 Auskunft. 


Tabellel 
| Doppelpunkt E Fußpunkte in der Nachbarschaft von E 

©, 1,672 333 2 1,672 333 3 1,672 334 1,672 4 1,6725 

I 3,396 159 = I 3,396 160 3,396 166 3,396 822 3,397 816 
mi 0.859 855 — 0,858 818 0,857 432 0,836 497 0,822 640 
m; Pe EEE 70.860.897 0,862 294 0,883 662 .. 0,898 287 
I ale 2,003 442 2,001 823 1,977 362 1,961 172 
I: 2.005 871 2,007 504 2,032 470 2,049 558 
ö1 12°4@/ 517° 12°46’8’ 12°45/117 12°30744’’ 1221.24 
Ör 7 2 12°47’34’’ 12°48’317 13342 13°12’48’ 
1/2(ö, + ör) 12°46’51” 12°46’51’ 12°46’51’ 12°46’52,5’’ 12°46’55’ 


In Tabelle 1 bezeichnen die Indizes Z und r an den Größen m, II, ö die Werte am linken 
bzw. rechten Fußpunkt der entsprechenden &, = const.-Kurve. J/ bedeutet den Quotienten des 
Druckes p, hinter dem Hauptstoß durch », vor ihm. Die letzte Zeile der Tabelle 1 zeigt, wie 
sich der Mittelwert 4 (ö,+ ö,) der Fußpunktwinkel ö, und ö, mit wachsenden ®, gemächlich 
nach größeren ö-Werten verschiebt. 


Der Doppelpunkt E gehört zu jener ‚Kurve‘ von konstantem @,, für welche die beiden 
Fußpunkte auf der ö-Achse zusammenfallen. Die Fußpunkte der Kurven befriedigen die 
Gl.(S)für2=0,d.h. A,= 0. Nunläßt A, die Faktorzerlegung 


4=(1+hf1+1— (1 — km] 


zu. Die beiden ersten Faktoren verschwinden im physikalisch sinnvollen Bereich nicht, so daß 
als Fußpunktgleichung eine Gleichung 4. Grads in m übrigbleibt: 


%=Pım'+P,m+ P,Rm+P,m+P,=0 ........(6) 
mit den Koeffizienten 
P,=1—-K% (+), PR,=—2(1— KR) k—(1—3%k)l+(1—k)1], 
P=—2(1—3)+4(1—k+2KR)1 + (—5+14k— KR —4R)® + (1 — KR, 
P,=2[2k+(1+k)l+2(—1+2%k— 2%)? +(1—3k)], R=(1+)[1—21+ (1—4k)l2]. 
Außer (6) befriedigt der Doppelpunkt E auch die Ableitung dieser Gleichung 99,/9m = 0, denn 
90 — 0 hat für 1 = I, die Doppelwurzel my. Nebenbei sei angemerkt, daß der Fußpunkt E’ der 


zum gleichen ®, = (®,)z gehörigen Kurve im ersten Lösungsbereich (s. Bild 2) bei öy, = 
— 21° 9 45’ liegt mit my — 3,084 045 und I/z, = 4,603 395. 


b) In Bild 4findet man in Kurve (, eine graphische Darstellung der Winkel öder Fußpunkte als Funktion 
der Anströmgeschwindigkeit o, für den 2. Lösungsbereich. Die Vergleichskurve (\, gibt den zu jedem w, ge- 


hörigen maximalen Ablenkungswinkel dar (©) im einfachen schiefen Verdichtungsstoß wieder. Für ihn . 


erreicht ö bei festem o, (bzw. !) ein Maximum, wenn m > 0 die Gleichung 
,=(l+k)mi (3 —I)m —-21=0. 2 Sr (7) 
befriedigt, wenn also 


_1—3+Y/®+2(1+4k)1+9 
— 2(1+k) 


‚ist. Setzt man dieses m in den Ausdruck für den Ablenkungswinkel im schiefen Stoß, z. B.in 


__U-krm(—m) ; 
(I+m[1+I1—(1-— km! u Su we © ee (8) 


m 
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ein, so erhält man Kurve (©, von Bild 4. Kurve (, begrenzt im 2. Lösungsbereichdas Gebiet der nicht- 
trivialen Lösungen des gegabelten Verdichtungsstoßes, d.h. sie begrenzt das Gebiet, in dem sich an die Lö- 
sungen für © =0 stetig solche mit 2 >0 anschließen. 
Kurve (, begrenzt das Gebiet der einfachen schiefen Ver- 
diehtungsstöße. 

c) Verfolgt man (z. B. in Bild 3) die sich im 
6, ö6,-Diagramm um E gruppierenden &, =const- 
Kurven für wachsendes o,, so findet man eine Kurve 
o, = (®,)r, welche die Gerade 6 — ö, in F berührt. 
Fist Doppelpunkt auf dieser Geraden, für die der 
hintere Nebenstoß zur Mac hschen Linie wird. Für 
sieist © —= m und die Gl. (5), die in diesem Sonder- 
fall vom 10. Grad in m wird, gestattet, den Faktor 
[1+2— (1 — k) m]? 5 [1+ (1+ k) m] vom 7. Grad 
in m herauszuziehen, so daß eine Gleichung 3. Gra- 
des für m übrig bleibt (I 18), nämlich 


RAM +QAm+QI=0 . . (9) 


mit 
,=(+ Mh}, 
9, =3—2k+3M%+2(—1+2k+ K)l, 
A=@—)U+k—(A—3M1, 


ee 
BEE 3 A Bild 4. Fußpunktkurve des 2. Lösungsbereichs C,. Maxi- 
Der Doppelpunkt F befriedigt außerdem die Ablei- male Ablenkung im einfachen schiefen Stoß 6, 


tung von Gl. (9) nach m. Beide Gleichungen sind 

simultane Bestimmungsgleichungen für die beiden Variablen =, und m = my des Dop- 
pelpunktes F. Successive Elimination von m führt zu 
IA 28-69 
2 Q no 9,9; Fi dd; ie Q,9; 
wobei das Gleichheitszeichen zwischen den beiden Ausdrücken eine Gleichung für I allein 
ergibt, da die Koeffizienten Q, nur von k und labhängen. Sie lautet für bel. k 


Dee. en u 1) 


My — 


e: (10), 


mit 
Ar = k?(1 —3k)?1%+Kk2(3 —18k+ 22 k? — 2 %K?+ 35 kt) 1° 
+k(3— 18 k+ 12 k? — 24 k?+ 122 kt-+ 72 k® — 27 k®) It 
+(1—6%k — 6 k2 — 60 k?+ 129 k?+ 290 k5 — 108 %®) 13 
— k (5 + 56 k — 5 k? — 438 k?+ 162 kt) 1? 
—2k(9+ 32 k — 147 kR?+ 54 kP)I—k (27 — 74k+ 272). 


Ar ist bis auf eine negative Konstante (— 4/27). k/[(1+ k)*(1 — k)®] die Diskriminante der 
kubischen Gl. (9). Für !< Ip ist Ay negativ und (9) hat nur eine reelle Wurzel m, die im 
1. Lösungsbereich den Endpunkt der betreffenden !—= bzw. ®, = const-Kurve auf der 45°- 
Linie x = m des Bildes 2 ergibt. Im ö, ö,-Diagramm müssen wir den Endpunkt und den 
Kreuzungspunkt einer solchen Kurve mit der ö —= ö,-Linie unterscheiden. Ein End- 
punkt entsteht in dem soeben betrachteten Sonderfall, ein Kreuzungspunkt im Fall ‚‚hinterer 
Nebenstoß ist senkrecht“. In Bild 2 dagegen liegen End- und Kreuzungspunkte des ö, ö,- 
Diagramms auf verschiedenen Kurven, die Endpunkte auf der Geraden & = m, die Kreuzungs- 
punkte auf der Kurve (',. Kreuzungspunkte gibt es nur im ersten Lösungsbereich. Den 
„Schleifen“ im ö, ö,-Diagramm (vgl. I Abb. 2), welche die Kurven des 1. Lösungsbereichs im 
Gebiet der gleichdrehenden Gabeln besitzen, entsprechen im z, m-Diagramm die Kurven- 
stücke zwischen C/, und der 45°-Linie «= m. Für 1> I» ist Ar positiv und (9) hat 3 reelle 
Wurzeln m. Für 1= 1, gilt (11) und 2 von den 3 reellen Wurzeln sind gleich. Für Luft ist 


Tabelle 2. 
Zu Bild 3 Zu Bild 2 en 
ö | 6 m ve: | Zu = IT, 
F} 0°55’16” 0) 0,053 673 0 1,062 711 a 
F 6°45°47° =öF 0,458 733 = my 1.535 983 
F, 26°24’12” 0 2,961 984 N) 4,460 780 | SE 1 
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Ir = 5,150 789, (@,)» — 1,814 881 und (10) liefert das zugehörige my (Tabelle 2). Die durch F 
laufende Berührungskurve habe F}, und F, als linken bzw. rechten Fußpunkt (s. Bilder 2 und 3). 
Tabelle 2 gibt über die Daten dieser Punkte Aufschluß. 

Die 3. reelle Wurzel von (9) für 1 = Ip ist my, = 3,775, öp = + 29° 13° 38”, sie gibt den 
Endpunkt F’ des dem 1. Lösungsbereich angehörenden Kurventeils auf der Geraden 2 —= m 
(Bild 2 und 6). 

Die in F berührende w, = (w,)p-Kurve ist bei wachsendem o, die letzte, die vom linken 
zum rechten Fußpunkt ohne Unterbrechung durchläuft. Die Bedingung x < m bewirkt, daß 
die &, = const.-Kurven des Bereichs 1,814 881<@,< 1,841 719 auf der Geraden = m 
(Bild 2) bzw. ö = ö„ (Bild 3) eine Unterbrechung er- 

,)- fahren. Die linken Fußpunkte dieser Kurven liegen im 
ni mo "225 &2 Intervall 0° 55 16”>6,> 0 (im Bild 4 auf Kurve (, 
zwischen F, und 0) und ihre rechten Fußpunkte im 
Intervall 26° 24° 12” < 6, < 27°43’ 16” (im Bild 4 auf 
Kurve (, zwischen F, und O,). In diesem Bereich füh- 
ren die 3reellen Wurzeln von (9) zu 3 Endpunkten der 
zugehörigen &, — const.-Kurve auf der Geraden = m 
bzw. ö = ö,. Einer davon gehört zum ersten Lösungs- 
bereich und liegt zwischen F’ und O, (Bild 2) auf der 
Geraden & = m. Die beiden andern gehören zu 2 Kur- 
venteilen des 2. Lösungsbereichs. Der linke Kurven- 
teil endigt auf der ö = ö,-Geraden in einem Punkte 
des Intervalls O0 <ö = ö,< ör (Tabelle 2). Aus Bild 5 
entnimmt man den zu einem bestimmten ®, gehören- 
den Endpunktals den kleineren der zweid = 6,- 
Winkel zur Ordinate &,.. Dergrößere Winkel zur 
gleichen Ordinate gibt das ö=ö, des Ausgangspunk- 
tes des rechten Kurventeils im 2. Lösungsbereich. Von 
dort läuft der rechte Kurventeil (Bild 3) ohne Unter- 
brechung über ein Maximum in ö,-Richtung bis zu 
seinem rechten Fußpunkt auf der ö-Achse. Die als 
Parameter an der Endpunktkurve (Bild 5) eingezeich- 
neten m-Werte deuten die zugehörigen & = m Werte 
der Endpunkte auf der 45°-Linie des 2. Lösungs- 
5? 1° 45m bereichs in Bild 2 an. 

Bild 5. Endpunktkurveim 2. Lösungsbereich zwischen d) Die obere Grenze des soeben betrachteten 
1,814 < o, < 1,842 @&,-Intervalls bildet eine ®, —= const.-Kurve, die den 

Ursprung 0 als isolierten Kurvenpunkt enthält (Ent- 

artung des linken Kurventeils). 0 liegt sowohl auf der m-Achse, als auch auf der & = m-Ge- 
raden (Bild 2). Dies bedeutet, daß sowohl (6) als auch (9) für 1 = I, die Wurzel m, = 0 besitzt. 
In beiden Gleichungen müssen also die Absolutglieder P, und Q, verschwinden. Sie enthalten 
in der Tat den gemeinsamen Teiler (1 —4%k)%?—21+1. Er hat=0 gesetzt, die beiden 
Wurzeln 7= 1/(1F 2yk). Diejenige mit dem Minuszeichen ist Z, daly>17>1 sein muß, — 
während die Lösung mit dem + Zeichen die in diesem Zusammenhang nichtinteressierende 
untere Grenze der / für Gabelstöße überhaupt ergibt. Für Luft wird I, = 5,578 157 bzw. 
(©), = 1,841 719. Für 1= I, reduziert sich die Fußpunktgleichung (6) wegen des Verschwin- 
dens von P,auf eine Gleichung 3. Grades für die übrigen Fußpunkte. Setzt man l,in g=0 ein, 


so läßt sich der Faktor (1 — Yk)? (1 -+ 2 Yk)? herausziehen und es verbleibt 
am? — 2 0, MEI Amer 8 Kb De 2) 


w; 


1,835 


1,830 


1,825 


1,820 


1815 
7,814 55 


mit 
“=(l+h(1+Yk)1—2Y), 8 =(1+4Ykfk, = —1+3%+10kyk+20 k2, 
Für Luft findet man die 3 Wurzeln m, = 3,437 080, m, —= 4,927 509, m; = — 0,936. Die erste 
Wurzel m, ergibt den rechten Fußpunkt 0, der (o,)» Kurve mit ö= 27° 43’ 16’ im 2. Lösungs- 
bereich. Die 2. Wurzel] m, führt zum Fußpunkt 0, im 1. Lösungsbereich bei ö = — 28° 41’ 36”. 
Die 3. Wurzel ist negativ, sie fällt daher außerhalb des Bereichs physikalisch sinnvoller Lösun- 
gen: O<m<l. 
Die Endpunkte der (o,)» Kurve auf der ö = ö,-Geraden folgen wegen Q, (l,) = 0 aus der 
quadratischen Gleichung Q,m?+ Q,m+ Q, = 0 für 1 = I, mit den Lösungen 
EN (13) 
2 (14h) (1+ Yh)* (1 —2YR) 
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worin D=1+5%k—12%k2+10(1—2%)yk bedeutet. Die Wurzel mit dem negativen Vor- 
zeichen gibt den Endpunkt 0, im 2. Lösungsbereich (Bild 2 und 3). Für Luft erhält man 
Mg, — 1,227 858, ö,, = (Öu)0, — 14° 24° 17, II, = 23,434 623. Die Wurzel mit dem positiven 
Vorzeichen gibt den Endpunkt 0, (Bild 2) im 1. Lösungsbereich. Für Luft wird m,, — 4,136180, 
ou = (d,)6 = 30% 8437, II, — 5,832 709 . 

e) Zwischen 0, und 0, liegt der Kreuzungs punkt 0, der (w,).-Kurve des 1. Lösungsbereichs mit der 


ö = öa-Geraden. Eine allgemeine Behandlung dieser Kreuzungspunkte findet manin I$ 5. Für unsere spezielle 
(@,)o-Kurve mit / =/, wird Gl, (I 21) 


4 
Bee ee senken « (14) 
j=0 
mit 
Re =(1+(1 +/R2(ı — 2/R%, RI = (1 +YhR(1 — 2/R) (3 — 8YE+k— 2ER), 
RD = _-X1 HH —2/H(2 +2Yk+3m, RO =-41+YH(ı—syH, RO <A. 
Für Luft folgen daraus die Daten für 0, (Bild 2), nämlich mo, = 4,693 690, ög, = (da)o, = 29°54’13, 
%, = 4, 008 398. 

Für einen Kreuzungspunkt gehört m zum starken schiefen Stoß, x zum schwachen schiefen Stoß für das- 
selbe. Daraus folgt C(x) = 0. Daß der hintere Nebenstoß trotzdem ein senkrechter ist, wird dadurch möglich, 
daß er sich nicht in der Anströmgeschwindigkeit w, befindet, sondern in einer kleineren, durch den vorderen 
Nebenstoß reduzierten 1» = 24. Aus der Bedingung des senkrechten hinteren Nebenstoßes folgt A(x) = 0. 
Beide zusammen A(x) = ((2) = 0 —s. $1(1) — führen zu Gl. (121). 

Beweis: Betrachten wir zwei m-Werte m, und m,, die bei konstantem / zu demselben ö im einfachen Stoß 
gehören, so gilt für sie nach (8) \ 

BERGEN mi (l — m,) . “| m3(l — m,) 
(1—-k2 [+7 —- (1 kKmP(l+m) M+I—(1-— km ]’(1-+m,) 
Diese Gleichung 3. Grads für eine der beiden Variabeln m, z.B. m,, reduziert sich, da sich der Faktor (1+ 1) 
(m, — m,) wegheben läßt, auf eine Gleichung 2. Grads 


m} bo (mı) — m, (1 — m,) [bo (m}) HL —K)mi] -Il+H)m(i—m)=0...... (1). 


Dabei ist b,(m,) die Größe b, aus (1) für m = m,. Setzt man m, = x, m, = m, so erkennt man, daß diese Glei- 
chung in C(x) = 0 übergeht, da c, = (l — m) [b, + (1 — k2) m?] ist. Aus (15) folgt sofort auch die Bedingung 
für jenes m, das zur maximalen Ablenkung der Strömung im einfachen schiefen Stoß bei festem I, d.h. o, gehört. 
Für dieses m ist m, = m, = m und (15) bat eine Doppelwurzel. In diesem Fall kann man aus (15) einen Faktor 
m, = m, = m herausziehen und es folgt 


mb —(— mb, +1 - Mm +1+1=0. 


Nach Einsetzen von b, aus (1) läßt sich diese Gleichung darstellen als das Produkt zweier Faktoren 
(1—-k)m — (1-+D][(l +k)m® +(3 — Im — 21] =0. 


Der erste Faktor ergibt eine Lösung m > I und fällt damit aus dem physikalisch brauchbaren Bereich heraus. 
Der zweite Faktor = 0 gesetzt, stellt die gesuchte Bedingung für ein zu Ömax(®,) gehöriges m dar und ist iden- 
tisch mit (7). j 

Für das oben betrachtete ®, = (w,), wird das zur maximalen Ablenkung im Hauptstoß gehörige m 


syE-ı +0 -VHV3 — Alk BE ri. (18). 
(1+%)(1—2/k) 

Für Luft ergibt sich Mpazx(ly) = 4,3844 und Ömax(@,0) = 30°20’57’. Es ist der höchste Winkel ö, ‚den die 

Schleife der (,)y-Kurve im ö, ö4-Diagramm zwischen 0, und 0, im Umkehrpunkt erreicht. Mmaz ist kein Maxi- 

mum der m-Werte, sondern m läuft monoton von 0,, wo esein Maximum m = I, hat, durch 0, über die Kurve Q, 

weg nach 0, auf der Kurve (w,), — siehe Bild 2. Für die Kurven des 1. Lösungsbereichs oberhalb (@,)z = 1,2448 

besteht, wie in $ 3 gezeigt wird, bei Luft kein Hindernis, in ihrer Schleife im gleichdrehenden Bereich Ömazx(®;,) 


zu erreichen, d, h, für sie ist maximale Ablenkung im Hauptstoß von gleichem Betrag wie im einfachen schiefe1 
Stoß möglich. 


f) Anders verhalten sich die Kurven des 2. Lösungsbereichs. Dies rührt daher, daß es 
zwei grundsätzlich verschiedene Ursachen gibt, weshalb eine solche Kurve im 6, ö.- Diagramm 
in ö-Richtung ein Maximum durchläuft. 1. die unter e) besprochene Ursache, nämlich daß 
öbzw.tg?öbeimonoton wachsendem m ein Maximum durchläuft. Dies ist ein Ausdruck 
für die Tatsache, daß es zu jeder im schiefen Stoß möglichen Ablenkung einen schwachen und 
einen starken Stoß gibt, die im Falle maximaler Ablenkung zusammenfallen. Wenn es die 
Gabelbedingungen zulassen, kann dieses Maximum der Ablenkung auch im Hauptstoß eines 
gegabelten Verdichtungsstoßes auftreten. Wenn sie es nicht zulassen, wird die Ablenkung im 
Hauptstoß jenes Maximum nicht erreichen können. Dann tritt 2. eine andere Ursache für ein 
Maximum, aber ein verringertes, auf, nämlich, daß mselbst, als Funktion von & (bei kon- 
stantem /) aufgefaßt, ein Extremum durchläuft, aber nicht den m-Wert erreicht, der 6 PR 
im einfachen schiefen Verdichtungsstoß gehört. Beispiele dafür geben für Luft die Kurven 
©, = 1,8 (maximales m — 2,8 im Hauptstoß einer Gabel, gegenüber m = 3,846 im einfachen 
schiefen Stoß), ©, = (w,)r, © = (@,), des 2. Lösungsbereichs (s. Bild 3). Bild 2 bringt den 


Mmaz (lo) = 
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Unterschied beider Ursachen klar zutage. Die Maxima der ersten Ursache liegen auf (,, die 
der zweiten auf C,. (Auch beim Überschreiten von 0, durchläuft m ein Maximum m = |, |6| 
durchläuft dabei ein Minimum, weil der z.B. in tg?ö auftretende Faktor I— m dabei ver- 
schwindet.) 

Für die Extrema der 2. Ursache gilt dm/dz = 0 neben Gl. (5). Aus (5) folgt durch Diffe- 
rentiation nach x unter Beachtung von m = m (%) bei konstantem !: 


6 6 
dm 9A 
n—1 Nn on — 
> n A,% er > an gr? 


n=1 n=0 


Für dm/d& = 0 erfordert dies 


6 
BE n A, gn—1 =E 0 ee EB Er LER 


n=1 


Bei vorgegebenem ! muß man % 
und mfür diese Extrema aus (5) 
und (17) als simultane Wurzeln 
berechnen. Für Fußpunkte die- 
ser Art, für welche nebenbei die 
zugehörige w,-Kurve ebenfalls 
senkrecht auf der ö-Achse des 
Bildes 3 steht, gilt außerdem 
& —0. Dann reduziert sich (5) 
auf, 4, = 0417) JauW AU 
Unter Benützung der Zerlegun- 
gen von A, und 4, bedeutet dies 
nach (115), = 0. und , =Ü: 
Ein solcher Fußpunkt tritt z.B. 
im 2. Lösungsbereich in der 
Nähe von w, = 1,75 auf. Aus 
diesem Punkt entspringt Kurve 
O,im Bild 2. 

g) Bild 6 enthält eine gra- 
phische Darstellung aller Kreu- 


45° 


J- d für Kurven I, V 
d. für Kurve 7 


(2) 
er 
2,0 2, 


0° 4 
WIE, 13 N 1ew Io 10 ATI AR LEER 2 120 er ©, —= const.-Kurven auf der 
Fr » ö&= Öö,-Geraden. Der Grenzfall 
Bild 6. Endpunktkurven 7 bzw. ZV im 1. bzw. 2. Lösungsbereich. Kreuzungspunkt- 1 ist 
kurve IT. 6,"Maxima im 2. Lösungsbereich auf Kurve /II. Zahlen in den Kiel „der hintere ‚Nebenstoß ist 
sen = Schnittzahlen der Herzkurven (s. $ 5b) für die vier o,, 6, Punktmengen links Machsche Linie‘ ist in Kurve I 
von I, zwischen I und IIT, zwischen III und IV, rechts von IV für den ersten, in Kurve IV für 


den zweitenLösungsbereich dar- 
gestellt. Kurve II gibt den Sonderfall des senkrechten hinteren Nebenstoßes wider. Für Luft 
gibt es demnach zwischen 1,1908 <o, < 1,814881 einen Endpunkt des ersten Lösungs- 
bereichs, zwischen 1,814 881 < o, < 1,8417 19 drei Endpunkte, einen des ersten und zwei des 
zweiten Lösungsbereichs und zwischen 1,841 719 <@, < @paz zwei Endpunkte, einen 
des ersten und einen des zweiten Lösungsbereiches für jedes w,, während es für 1,353 37 <o, 
< Omas Immer einen Kreuzungspunkt im ersten Lösungsbereich auf der ö = ö,-Geraden 
gibt. Die Kurve III gibt den maximalen ö,-Winkel des 2. Lösungsbereichs zwischen 
(©,)z — 1,6723 < @, < (@,)x — 1,9936 wider, worauf in $5b näher eingegangen werden 
wird [vgl. Linie EX in Bild 3]. 


$ 3. Der erste Lösungsbereich unterhalb o, = Y2 — k 

Außer der Kurvenschar des 2. Lösungsbereichs gehört bekanntlich auch ein Teil des 
1. Lösungsbereichs zu nur gleichdrehenden Gabeln. Sie gruppieren sich in der Nähe des Null- 
punkts im ö, ö,-Diagramm. Ihr I bzw. w® liegt in den Grenzen 1/(1+2 Yk) <I<lY(1—k) 
bzw.2/1+Yk)<o?<2—k. Für Luft ist dies der Bereich 1,190 83<w,< 1,353 37. Die 
untere Grenze gehört zum tiefsten &, für Gabelstöße überhaupt, die obere zum tiefsten o, für 
eh mit senkrechtem Hauptstoß. In Bild 7 ist eine Folge solcher &, = const.-Kurven 
skizziert. 

Kurve II in Bild 8 zeigt für Luft, wie der Winkel ö der Fußpunkte der ,-Kurven dieses 
Bereichs mit wachsendem o, von ö = O0 für @, = 1,1908 bis zu einem Maximalwinkel im Punkte 
M anwächst. In M ist 1= 0,771 188, (w,)»r = 1,252 041, m = 0,594 223, öyr = 7° 20 467, 


zungs- und Endpunkte von: 
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II = 1,694 289. Aber ö,, ist nicht zugleich ein maximaler Ablenkungswinkel im einfachen 
Verdichtungsstoß für &, = (®,)y , den wir Öyaz (@,,,) nennen. Denn aus (7) folgt für obiges o, 
der Wert m = 0,539 459, woraus Öyaz (@,,) = 7° 29 50” hervorgeht. Da mals Funktion von & 
entlang der ®, = (®,)y-Kurve von einem größeren m = 0,594 für 2=0 monoton nach 
einem kleineren m = 0,168 für 
x —= m (45°-Linie) läuft, nimmt es 
auf der Kurve auch den Wert, der 
ZU Öyas gehört, an. Der maximale 
Ablenkungswinkel wird also in 
Übereinstimmung mit $2f. im 
Hauptstoß erreichbar. Für &, = 
(©,)ır geschieht dies aber nicht auf 
der ö-Achse (ö, = 0), sondern in 
einem Punkt ö6,<0 des ö, ö,-Dia- 
gramms. Die ©, = (o)y-Kurve 
steht also auf der ö-Achse nicht 
senkrecht, sondern verläßt sie mit 
wachsendem ö, in positiver ö-Rich- 
tung. 


Für ®,> (o,)y nimmt der 
Fußpunktwinkel ö wieder ab, um 
bei der oberen Grenze des Bereichs 
©, = 1,3534 den Nullpunkt noch- 
mals zu erreichen und für noch 
höhere Werte von o, endgültig zu 
negativen ö- Werten überzugehen. Bild 7. Kurven des 1. Lösungsbereichs für kleine w, 

Zu jedem Fußpunktwinkel im Be- 

reich O<ö< (ö)y, gehören zwei w-Kurven dieses Bereichs, eine mit kleinerem und eine mit 
größerem @, als (®,)yr. Die beiden Werte entnimmt man aus Bild 8, wo sie zur gleichen Ordi- 
nate der Kurve II gehören. Nur in öy selbst mündet die Kurve ®, = (®,)y allein. 


Dabei ist die in $2 behandelte Kurvenschar außer acht gelassen. Zieht man sie noch in 
Betracht, so vermehrt sich die Zahl der Kurven, die in einem Fußpunkt münden, um eine. 
Zwischen O<ö< (ö)y haben wir das Achsenstück des Bereichs der 3-fachen Überdeckung des 
ö, ö„-Diagramms (s. Fr. Wecken [6] Bild 17) vor uns. 


Für welches &, steht 
die &,-Kurve senkrecht auf 
der ö-Achse? Darüber gibt 
die Kurve III des Bildes 8 
Aufschluß. In ihr wurde der 
maximale Ablenkungswinkel 
für den einfachen Stoß gra- 
phisch dargestellt. Sie be- 
rührt die Fußpunktkurve II 
im Punkte B, dessen !- und 
m-Werte also simultan die 
GlIn. (6) und (7) befriedigen. 
In B ist 1 = 0,743 428, (o,)» 
— 1,244 768, m = 0,519286 
und 05 = 7° 9.14”. "Für 
©, = (®,)z wird der maxi- 
male Ablenkungswinkel ge- 
rade auf der ö-Achse ange- , | 
nommen (B in Bild 7). Die 14 1.20 1.22 724 126 128 1,30 1,32 1,34 1,36 (©), 
©, = (®,)g-Kurve steht hier Bilds. Endpunktkurve 7, EN II, ER III für kleine o, 
auf der ö-Achse senkrecht. ir. 
öz ist 11’ 32” kleiner als öy. Für &,> (®,)z rückt der maximale Ablenkungswinkel öy.. (@,) 
auf der gestrichelten Linie in Bild 7 aus B hinaus zu höheren ö,, während er für ®, < (@,), im 
Hauptstoß eines gegabelten Verdichtungsstoßes nicht erreicht werden kann. Für ®, zwi- 
schen (9,)s < ®, < ®mas wird jedoch der Richtung des Hauptstoßes im 1. Lösungsbereich 
keine stärkere Einschränkung als im einfachen schiefen Stoß auferlegt. 
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Die &,-Kurven für &,<(@,)z laufen unter einem stumpfen Winkel, diejenigen für 
©, > (w,), unter einem spitzen Winkel aus der ö-Achse im Öö, ö,-Diagramm aus (Bild 7). 

In Kurve I des Bildes 8 sind schließlich die ö — ö,- Winkel der Endpunkte der Lösungs- 
kurven zwischen 1,1908 < &, < 1,3534 auf der ö = ö,-Geraden aufgetragen. Der Kreuzungs- 
punkt der Kurven I und Ilin Bild 8 bedeutet z. B., daß die ®,-Kurve für w, — 1,289 von der 
ö-Achse bei 6 — 5° 28’ ausgeht und beim gleichen Winkel 6 = ö, = 5° 28’ auf der ö = ö,-Ge- 
raden endigt. Der zu ®, = 1,289 auf Kurve III gehörende Wert von Öyas (1,289) — 8° 58° 
gibt das höchste dan, das zwischen Fuß- und Endpunkt erreicht wird. Die Kurven des Bildes 7 
beginnen bei ihrem Fußpunkt auf Kurve Ilin Bild 8 und endigen bei ihrem Endpunktwinkel 
auf Kurve I des Bildes 8. Aber die Kurven für ®, > (o,)z laufen von Kurve Ilerst nach Öyar 
auf Kurve III und kehren von dort nach I zurück. Jedes ö der Fläche, die in Bild 8 von I, III 
und dem rechten Teil von II zwischen I und III begrenzt ist, wird zweimal durchlaufen. 

Die aus dem Punkt M in Bild 7 auslaufende gestrichelte Kurve ist die Enveloppe der 
Kurvenschar ®, < (@®,)yr. Ineinem Punkt dieser Art hat ö bzw. tg? ö[Gl. (8)] als Funktion von 
lund m bei festem ö, ein Maximum (dö—0) und ö, als Funktion von Z und & bei 
festem ö ein Minimum (dö,„= 0). Für den Zuwachs dö des Ablenkungswinkels ö als Funk- 
tion der Zuwächse dl und dm der Variablen l und m folgt aus (8) 


(4 —k){m(1+ m)[1 —1+(1+ k) m]d—(1+ I)[( + k) m®—(l— 3) m — 2 1]dm} 
e: 2(1+Db,Yl — m) (1-+ m) 

Für ein Extremum ist dö—=0. Dafür muß der Zähler von dö verschwinden, da der Nenner 
für endliche 1 > m> O endlich bleibt. Dies ergibt für ein solches Extremum die Bedingung 

dl 1+1 (1+%k) m? — (l—3) m — 21 
dm 1-+m m [1 —1+ (1+ k) m] 


Da M auf der ö-Achse liegt, müssen seine l, m-Werte neben (19) auch Gl. (6) 9, = 0 befriedigen. 
Durch Differentiation von 9, = 0 nach m folgt 


au = dr, m" Nap mri—h( (20) 
dm + a n ZU BE SEN A ER TE REERER K 


ı=0 


dö 


(18). 


HEN. 


Ersetzt man darin dl/dm durch die rechte Seite von (19), so entsteht eine weitere algebraische 
Gleichung, die durch die Zund m-Werte von M simultan mit 9, = 0 erfüllt werden muß. 

Für einen Punkt der Kurve, die in Bild 7 aus M als gestrichelte Linie ausläuft, gelten 
1. die Gl. (5) und 2. die daraus durch Differentiation nach !, m, & bei konstantgehaltenem 6, 
entstehende Gleichung 


6 6 
0A, , 9A, dm\ ,, da » 
N — A,ar 1 —= 6. Green Eee e 
>| 1 "m zE ee a go 
Für dm/dl ist der reziproke Wert der rechten Seite von (19) in (21) einzusetzen. dax/dl findet 
man aus der Nebenbedingung dö,—=0. Da tg? ö, von x und Lin gleicher Weise abhängt wie 
tg? ö von m und ], so erhält man da/dl für dö, = 0 aus dm/dl für dö — 0, indem man in (19) m 
durch z ersetzt. Führt man also noch 


02° 1% [11 —!+(1-+Kkal 
u 17 ed 
in (21) ein, so hat man in (5) und (21) zwei simultane algebraische Gleichungen für I, m, z, 
welchen die Punkte der aus M in Bild 7 auslaufenden gestrichelten Kurve genügen müssen. 
(21) reduziert sich für x = 0 wegen da/dl = 0 auf (20) für den Kurvenpunkt M. 
Wird dagegen I konstant gehalten, so ist dl/dm = 0, und (19) geht in A, = 0 Gl. (7) über. 
Punkte auf der gestrichelten Kurve, die in Bild 7 aus B ausläuft, befriedigen (5) und (7). 


$ 4. Der Grenzfall maximaler Anströmgesehwindigkeit für Luft 

Wenn die Anströmgeschwindigkeit die höchstmögliche &, = @yar ist, muß man 2 Fälle 
unterscheiden, die in1$6 als Falla) und Fall b) gekennzeichnet wurden. Im Fall a) ähneln die 
Strömungsbilder denen für geringere Anströmgeschwindigkeit. Ihr Diagramm wird jedoch 
dadurch etwas vereinfacht, daß sich (5) vom 6. auf den 4. Grad reduziert — siehe Gl. (125). 
Im Fall b) tritt etwas Neues, nur im Grenzfall verschwindenden Druckes in der Anströmung 
Mögliches auf. Der vordere Nebenstoß nimmt nämlich die Richtung der Machschen Linie an, 
ohne daß das Druckverhältnis J/, = 1 zu sein braucht. (siehe Fr. Wecken, Bild 10 und 11). 
Diese Mac hsche Linie unterteilt das Feld der Anströmung in zwei Felder z. B. vom Druck 0, 
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deren Druckverhältnis 0/0 jedoch im allgemeinen von 1 verschieden ist. In diesem Fall redu- 
ziert sich (5) auf das Quadrat eines in «linearen Ausdrucks. s. Gl. (129). 

Als Anwendung unserer Gleichungen 125 und I 29 berechnen wir für Luft die zuo, = Oyar 
gehörende Kurve im ö, ö,-Diagramm in 34 Punkten. Das Ergebnis ist in Bild 9 verwertet. Die 
Punkte 1 bis 7, 34 bis 28 liegen auf 2 Strecken der ö-Achse symmetrisch zu 6 —=(, sie sind 


ah, DT ae] 


40° 


11234 5 6, 
Be 
RE ee eo ae 


Bild 9. ö, ö,-Diagramm für o, = Omar 
Lösungen für Fallb. Die Punkte 7 bis 19, 20 bis 28 liegen auf 2 Kurventeilen als Lösungen des 
Falles a. An den Punkten 7 und 28 hängen beide Fälle zusammen. 


Für &, = Oyaz werden I, m— mw. Im Falle a) wird außerdem auch 2 —= ®. Aber in den 
Grenzwerten ihrer Quotienten « = lim m/l und &= lim «/l gilt obenerwähnte Gleichung 


l, mo» l,x> © 


4 
4. Gradsin &E mit den von u abhängigen Koeffizienten B,, & B„&"—=0 (125). 

: n=0 

Im Falle b) bleibt xjedoch mit Ausnahme des Punktes 7 endlich. Im Punkt 1 beginnt er mit 

z —= 0 und u (0) = 1/(1-+ k) = 0,8559. Unter Durchlaufen der Punkte 2 bis 7 wächst x und 
erreicht in 7 den Wert &, wobei u—= u(z) Gl. (131) langsam auf « (o) = 0,9943 zunimmt. 
Dasselbe gilt symmetrisch zu ö = 0 zwischen den Punkten 34 und 28. Längs beider Geraden- 
Braske autder Achse vonsBild. Ist 75 055, = 0 "N =8, In =, md), 
6,—=6. Über die übrigen Größen gibt Tabelle 3 Aufschluß. 


Tabelle 3 

1/34 | 0,8559 0 45°22/37.° 1 

2/33 | 0,9 0,1602 44°45’44.’ 1,1871 
3/32 | 0,9334 0,3904 42°48’46’’ 1,4561 
4/31 0,95 0,6127 40°47’57’° 1,7159 
5/30 | 0,975 1,6100 SEIT BTE 2,8812 
6/29 | 0,99 7,5825 95595. 9,8594 
7/28 | 0,9943 oo 19°50’ 9” (6) 


Das obere bzw. untere Vorzeichen vor ö in Tabelle 3 bezieht sich auf die erste bzw. zweite 
Punktfolge in der ersten Spalte. Der vordere Nebenstoß ist in diesem Falle zugleich Strömungs-, 
Trennungs-, und Machsche Linie. Die Formel für die in I$ 1 eingehender definierten Größen 
vereinfachen sich hier unter Benützung des Grenzwertes u zu 


_B-A-hulll+hu—1] > en 
AI (l+hutrkkli— Bd) 
lim. © H A DL: a: 


animal Itlitke’ viel (A+ka)li+li+Ne 


u(1+ kx) 
tg? 9,, — ketg? 9,,. 
1+2— (1+ka) u’ Bo ur 
Dabei bedeuten #’, U’ die zu x, lJanalogen Größen im Strömungsfeld vor dem hinteren Neben- 


stoß. In den Bildern 10 bis 12 findet man eine graphische Darstellung der Winkel zwischen den 
Strömungsrichtungen und den Verdichtungsstößen, 


tg?9, = er ;, t9,—=kttg?d; td, = 
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Fall a) ist mannigfaltiger. Auch er zerfällt in 2 Teilstücke im 6, ö,-Diagramm. Das erste 
Stück erstreckt sich von Punkt 7 bis 19, das zweite von 20 bis 28in Bild 9. Für die Endpunkte 


129° 


| een 


50° n er : 
50 EN 30 20 O0 E20 E30 a 0° TOE EEE E FOIM SO 
Bild 10. Winkel d, zwischen Anströmung und Hauptstoß für Bild 11. Winkel da zwischen Anströmung und vorderem 
©, = Opa: Pı läuft von Punkt 1 bis 19, springt nach 20, läuft Nebenstoß für ®, = ® 


von 20 zurück bis 28, kehrt dort um und läuft bis 34 _ 


\ 7) i 
280 60 30% 2070 WO 102 EOS 03 
—/ a 


Bild 12. Winkel 9,5 zwischen Strömung vor dem hinteren Nebenstoß und 


dem hinteren Nebenstoß für o, = Omar 


19 und ist &=u. Für 19 ist m =&= 1/(1+ k) = 0,8559, 6 = &, — 45° 22° 37” und für 
20 ist = E= (1 — 3 k)/(1 — k)? = 0,7155, 6 = &, = 42° 48° 46” für Luft (siehe 127). Die 


Tabelle 4 
Bet £ ö du I; 

8 0,997 0,0102 —14°54’10”’ 4°49'29’ 97,441 

9 0,999 0,0474 — 8°49’43’’ 10°25’15’’ 21,087 
10 1 0,1386 0° 17°59’43’’ 7,213 
11 0,999 0,2608 8°49’43’’ 24°59'55 3,830 
12 0,997 0,3560 14°54’10’ 29°29’22” 2,801 
13 0,9943 0,4405 19°50’ 9° 33° 572" 2,257 
14 0,99 0,5251 25° 5’25” 36°23745 1,885 
15 0,975 0,6800 34°27737’ 41°45’28’ 1,434 
16 0,95 0,7794 40°47’57’ 44°25’21’’ 1,219 
17 0,9334 0,8101 42°48’46’’ 44°59’12° 1,152 
18 0,9 0,8411 44°45’44 45°19°47’’ 1,070 
19 0,8559 0,8559 45°22'37° 45°22’37’ 1 
20 0,7155 0,7155 42°48’467 42°48’46’’ 1 
21 0,7967 0,6396 44°45’44’ 40°27’427 1,246 
22 0,8559 0,5557 45°22'37’’ 37°31/51” 1,540 
23 0,9 0,4594 44°45’44’ 33°50’48” 1,959 
24 0,9334 0,3480 42°48’46” 29° 8713’ 2,682 
25 0,95 0,2716 40°47'57” . .26°32’13” 3,498 
26 0,975 0,1165 34°27’37' « 16°27’16’ 8,373 
27 0,99 0,0132 25° 5/25’ 6°28’31” 75,175 
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Winkeldifferenz beider Endpunkte ist A6 = A 6, = 2° 33’ 51”. Auf sie wurde schon von 
früheren Autoren aufmerksam gemacht. Die Daten der Punkte 8 bis 28, für die durchweg 
x—=m,/7=II, = xgilt, wurden aus der Gl. 4. Grads für £ (125) berechnet und in Tabelle 4 
zusammengestellt. 


Die Größen x’ und /’im Strömungsfeld vor dem hinteren Nebenstoß bleiben im Fall a) 
endlich: 


a te & 
 (1+k)E 

Ba (HE = 
k(+k)E 


und für die Winkel zwischen 
den Strömungsrichtungen und 
Verdichtungsstößen gilt 


/b) 
u > b 
£ 
tg? At ri weg F 
1+r 
tg? 020) Zr Bez 
1—E—ku' 


Korn kit’, 
Pe 

ku\% 
tg? d 6-2) tg? 9,2 
Betrachten wir noch im Bild 9 
Kurvenpunkte, die zur gleichen 
Abszisse, also zum gleichen Ab- 
lenkungswinkel im Hauptstoß 
gehören. Einige unter ihnen 
stehen in der Beziehung (15) 
zueinander. Für 0, = max 
d.h. 2= oo, vereinfacht sich 
diese Gleichung, wenn man sie (6) (h) 
mit /* dividiert und zur Grenze 
m,l—= oo übergeht. Es sei 


Punkt 26 


(8) 


im m 

a = lm’ —, = lim — 

m, ,‚I>® l 

dann wrd lim ———= 

BER U 

[1 —(1— K?)u,]. Das Absolut- 

glied geht gegen 0, Gl. (15) 

reduziert sich auf eine lineare, 
deren Lösung 


1. 
BT —(1 Ba} kı (26) Punkt 30 
lautet. w, und zw, stehen zu- 
einander in der Beziehung von 
schwachem und starkem schie- 
fen Stoß zum gleichen Ablen- 
* kungswinkel. Diese Beziehung (b) 
besteht zwischen den Kurven- Bild 13, Strömungsbilder für die Punkte 15, 26, 30 des Diagramms Bild 9 
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stücken 20 bis 22 (schwacher) und 17 bis 19 (starker Hauptstoß). An den Punkten 19 und 22 
wird Öyaz (Omas) erreicht, dort ist 45 = u = Y(1+k). Beziehung (26) zwischen den beiden 
Punkten 20 und 17 ergibt den u-Wert für Punkt 17, da derjenige von Punkt 20, u, = 
(1— 3 k)/(1—k)?, bekannt ist. Man erhält u,—= (1+ k)/[(1—k) (1+3k)]für Punkt 17. Punkte 
außerhalb dieser Strecken, die zum gleichen ö, aber zu verschiedenen Kurvenstücken des 
Falles a) oder b) gehören, also Punkte zwischen 13 und 19 verglichen mit solchen zwischen 28 
und 22 oder solchen zwischen 28 und 34 haben gleiches u für gleiches ö. Drei Punkte dieser Art 
sind 15, 26, 30. In Bild 13 findet man die Strömungsbilder für diese Punkte, sie haben den- 
selben Ablenkungswinkel im Hauptstoß und gleiche 9. 


$ 5. Algebraische Behandlung der Herzkurvenmethode 


a) Die Herzkurvenmethode von A. Weise [2] bringt zur ne der beiden Gabel- 
bedingungen ö = 6,+ 6,, /T= II, : II, in einem 6, II-Diagramm 2 Herzkurven zum Schnitt, 
die Herzkurve für den Hauptstoß (Hauptherz) mit derjenigen für den hinteren Nebenstoß 
(Nebenherz). Die Gleichung für das Hauptherz erhält man aus der Gleichung für den Ab- 
lenkungswinkel ö als Funktion von of, o,, und k 


© 


= nn eat 2 
(el Kae 


On, 1+(1—k)ol —o,, 


indem man darin statt ©} und @,, die Druckverhältnisse //, und I/ einführt. 


_ ih) 
’ (1 kob) 


bedeutet das Druckverhältnis (p,/p,),s im senkrechten Hauptstoß, worin das Druckverhältnis 
im schrägen Hauptstoß // = p,/pı = (wo, — k)/1— kw;,) für ©, = w, übergehen würde. 
Damit nimmt die Gleichung des Hauptherzes die einfache Form an 


= Made) Henn 
1 (DIRT NDEE 


II, ist der höchste Druck im Hauptherz, d. h. die Koordinate seines Scheitelpunktes, während 
für seine Spitze /7’=1 gilt. Um den Nullpunkt des Koordinatensystems in die Spitze des 
Herzes zu verlegen, die Kurven in der Ordinatenrichtung gedrängter zu gestalten und die so- 
fortige Ablesung des Druckverhältnisses //, im Schnittpunkt von Haupt- und Nebenherz zu 
ermöglichen, wird man in der zeichnerischen Praxis wie Wuest ([5] Bild 2) auf der Ordinate 
logarithmischen Maßstab A = log // usw. benützen. 

Die Herzkurvenmethode greift nun ein bestimmtes Druckverhältnis //, im vorderen 
Nebenstoß heraus. Dazu gehört ein bestimmter positiver Ablenkungswinkel ö,, den man auf 
dem Hauptherz ablesen kann. Analytisch bedeutet dies die Berechnung des Winkels ö, aus 
127) Aurel: 

ji ig ö Ye AUF) IL 1]; 
* 1—-(1—kI,+I,Y I,+k 


In dem durch //,, ö, markierten Punkt des Hauptherzes legt man die Spilze des Nebenherzes 
mit gleichgerichteter, paralleler /7,;-Achse. Die Bestimmung seiner Form erfordert die Kenntnis 
der Anströmgeschwindigkeit ® = 4 = @,, des hinteren Nebenstoßes. Dies bedeutet in 
unserer Darstellung die Kenntnis des Druckverhältnisses IT,, für den senkrechten hinteren 
Nebenstoß 


[0 (47): 


WAR 


__..@1b k a 1 (II, + k) (I, — II,)\ 
ne lH kl; 


Diese Formel erhält man, wenn man in 


. (29). 


Y— On, , 1 1—kat 

1— kon, mul —kon. 

die Druckverhältnisse /7, und /Z, — (on, „ — k)/1—k On.) einführt. Dann wird in der ersten 
Gl. (29) 


3 2 
0 = Wu = 


a + k)UR—I)+ (LH EIR) 
A (1+ kIL,) (I, -+ k) 
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Mit I7,, lautet nun die Gleichung des Nebenherzes (Gl. (27) mit IT, statt IT, IT,, statt IT,): 


(1 =) Vs) II,; —I, 
1—(1— kA), +2,fY I,+k 
Wenn man Zähler und Nenner mit IT}? erweitert und die zweite Gabelbedingung 
IT—= II, : II, beachtet, wird daraus 
DE SIT IT, TEE TI else = wel in 


. (30). 


29%, = 


tg Ö5 = 


Das Produkt I/,, IT, ist die geschweifte Klammer in (29). Setzt man sie statt /T,, IT, in tg 6, 
ein, so werden alle auf den hinteren Nebenstoß bezüglichen Größen eliminiert. 


Die Anzahl der Schnittpunkte von Haupt- und Nebenherz hängt bei vorgegebenem k und 
II,, d.h. bei vorgegebenem x und vorgegebener Anströmgeschwindigkeit @,, von der Wahl von 
II, ab. Die Unbekannte, die durch einen Schnittpunkt beider Herzkurven bestimmt wird, ist 
II, das Druckverhältnis im Hauptstoß. Damit ergibt sich auch ö und 7, = II/II,, 6, = 6 — Ö,. 
Analytisch bedeutet das Aufsuchen von Schnittpunkten die Anwendung der ersten Gabel- 
bedingung ö=ö, + ö,. Sie führt unter Beachtung von (27), (28) und (31) zu der Gleichung 


[ AAN) REN _ ger AZBUR—N /IR—IL) 

ee erh Eat, -A—naN Inte) 
[|  a—hW(M—IL) 1/0. — I 

ee en Treo 


arcig 
(32), 


worin noch //,, Z/„ durch die geschweifte Klammer in (29) zu ersetzen ist. Mit Hilfe des 
Additionsgesetzes der arctg und 2-maligem Quadrieren folgt daraus eine algebraische Gleichung 
zwischen %, /7,, II, II, mit rationalen Koeffizienten. Ohne sie explizite aufzustellen, kann 
man über sie aussagen, daß der für den gegabelten Stoß wesentliche Bestandteil in ihr vom 
6. Gradin jeder der Variablen /7,, IT, IT, sein muß. Denn diese Variablen /,=1-+(1+k)l, 
I=1+(1+k)m, I,=1+(1+k)x hängen linear mit den Variablen !, m, x unserer 
Gl. (5) zusammen und diese treten dort alle gleichermaßen in höchstens 6. Potenz auf. 


b) Wir wollen für Luft die Anzahl von Schnittpunkten von Haupt- und Nebenherz angeben, die für ein 
bestimmtes w, auftreten können. Da beim Anlegen des Nebenherzes I/, vorgegeben wird und damit bei vorge- 
gebenem o, auch ö„ bestimmt ist, bedeutet dies die Bestimmung der Anzahl von Schnittpunkten, die im ö, 
öa-Diagramm (Bild 3) eine Gerade ö, —= const., also eine Parallele zur ö-Achse in der Höhe ö,, mit einer 
©, = const.-Kurve besitzt. 

Eine solche Parallele hat einen Schnittpunkt mit einer bestimmten &,-Kurve des 1. Lösungsbereichs, so- 
lange ö, kleiner ist als das ö, des Endpunkts dieser Kurve auf der ö = ö,-Geraden des Diagramms. In Bild 6 
ist in Idie Endpunktkurve des 1. Lösungsbereichs wiedergegeben. Links von Kurve / hat manalso keinen, 
rechts davoneinen Schnittpunkt im ersten Lösungsbereich. 

Im 2. Lösungsbereich gibt es zwischen (w,)z < ©, <(o,)r zwei weitere Schnittpunkte der beiden 
Herzen, wenn das ö, der Parallelen zur ö-Achse kleiner als das maximale ö, der w, = const.-Kurve ist, und 
keinen, wenn es größer ist. Diese maximalen ö, wurden in Bild 6 als Kurve I/II zwischen E und X eingezeichnet, 
als Kurve C, in Bild 2, als Linie EX in Bild 3. Im Intervall (o,)r < ®, < (o,), erfährt die w,-Kurve eine Unter- 
brechung auf der ö = ö,4-Geraden. Wenn die Parallele diese Unterbrechung durchschneidet, gibt es statt 2 nur 
einen weiteren Schnittpunkt der Herzen. Liegt die Parallele tiefer als der untere oder höher als der obere End- 
punkt auf der ö = ö,-Geraden, aber tiefer als das Maximum der w,-Kurve, so gibt es wie früher 2 weitere Schnitt- 
punkte, Die Kurven /V bzw. III in Bild 6 begrenzen also das Gebiet, in dem ein bzw. zwei weitere Schnitt- 
punkte der Herzen auftreten. Oberhalb (w,), existiert kein unterer Kurventeil im ö, ö„-Diagramm mehr. Dort- 
tritt ein weiterer Schnittpunkt auf, solange die Parallele unterhalb des Endpunkts der &,-Kurve auf der ö, ö,- 
Geraden liegt; zwei weitere Schnittpunkte, wenn sie zwischen diesem Endpunkt und dem Maximum der Kurve 
durchläuft; und kein weiterer Schnittpunkt, wenn sie oberhalb des Maximums liegt. &, = (w,)x = 1,9936 ist 
das höchste w,, bei welchem ein solches Maximum auftritt, esist zugleich Endpunkt der Kurve auf der ö, ö„-Ge- 
raden für den 2. Lösungsbereich. Oberhalb (o,)x gibt es also entweder keinen oder einen weiteren Schnittpunkt 
der Herzen im 2. Lösungsbereich, wenn die Parallele höher oder tiefer als der Endpunkt der untersuchten o,- 
Kurve auf der ö, ö„-Geraden liegt. 

Bild 6 gibt zusammenfassend über die Anzahl der Schnitte von Haupt- und Nebenherz Auskunft. Man 
greife einen bestimmten Punkt w,, ö„in Bild 6heraus. Liegt erlinks von Kurve I, sogibteskeinen Schnitt- 
punkt; liegt er zwischen den Kurven / und III bzw. zwischen I und IV oberhalb Punkt X, so gibteseinen 
Schnittpunkt ; liegt er zwischen den Kurven /// und IV, so gibt es drei Schnittpunkte und liegt er rechts von 
Kurve IV, so gibt esz wei Schnittpunkte. Die durch Kreise umrahmten Zahlen in Bild 6 deuten diesen Sach- 
verhalt an. Kurve I/ spielt bei dieser Betrachtung keine Rolle. Diese Ergebnisse stimmen mit denen von 
Fr. Wecken überein. 


c) Wir berechnen noch die Daten des Punktes X, d.h. die obere Grenze der Dreischnittig- 
keit der Herzkurven. X liegt einerseits auf der Geraden & = m, andererseits auf der Kurve Q, 
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der Maxima in Bild 2, wofür da/dm = 0 gilt. Differentiert man (5) nach m unter Beachtung 
von = x (m), so erhält man mit da/dm =O0 und = m 


= om 


Diese Gleichung ist vom 9. Grad in m; sie läßt jedoch den im physikalisch interessierenden 
Bereich nicht verschwindenden Faktor 2 [1+ 1— (1— k) m] b, vom 3. Grad in m herausziehen. 
Dann bleibt eine Gleichung 6. Grads in m übrig, die man simultan mit Gl. (9) für die End- 
punkte auf der Geraden x = m erfüllen muß. Durch Kombination beider Gleichungen läßt 
sich die erste auf den 4. Grad in m reduzieren, ohne daß dabei der in %k, } ganze, rationale Cha- 
rakter der Koeffizienten verloren geht. Man zieht von der ersten Gleichung 


{1 —k2(1+k)2+3k)m’+ [2 —10k—5%K2+9%k) + (—3 —3%k+ 2%) 1]m?} 
mal Gl. (9) ab und erhält dabei folgende Gleichung 4. Grads in m: 


4 
2 B,m 0 N re 


n=0 
mit den Koeffizienten 
H, = —3+4k+27k? —3k? — 18 k?+9%k5+(—4+8k+22K2 +23 +2%°+2%5)l 
+1+8%+7 RB +43), 


H, = —7+15k+41k2 — 49 k®+ 30 kt + (—5+ 22 k — 23 k?+ 32 k®-+ 32 kt) ] 
+(12 —21k— 172 +31 R3+ 17 kK?+AKkd)P+(— 2 +2k+ Kr —ARB + Kt), 


H,— —8—k— 9k?+ 36 k®+ (1—19 k—54 k?+ 108 %3) + (18— 50 k+ 15 k?+21k°+ 12 ke) 12 
+ 10420 ® +k— 17 + 2K)B+(L—2K— B)l, 


H=—5—3k+18hK+ (4— 106432 2) 14 (7 —4k+ 39 k? — 24 K8) 12 
+(—104+Bk—6R—UR)B+A@—5k— Me), 


HH DAS Bi gel ame 


Aus dieser Gleichung zusammen mit Gl. (9) sind die Daten /y und my des Punktes X zu 
berechnen. Für Luft erhält man !y = 8,9948, my —= 4,100, (o,)x = 1,9936, (@,,)x — 1,7369, 
6 —= 6, = Ör —= 26° 53’ 45°, IIx = 5,7904. Die Kurve, die in Bild 3 die Maxima zwischen Z 
und X verbindet, kann man in guter Näherung als Gerade ansehen 

Öx 


6, & —— (6 — ö5) = 1,9055 (6 — 65) , 
Öxr — ÖR 


die den Winkel 62° 18,5’ mit der ö-Achse bildet. 


Der eine von uns (Th. Vollmer) ist Herrn Professor Dr. A. Weise für die Über- 
tragung einer Diplomarbeit auf diesem Gebiet zu großem Dank verpflichtet. Ferner danken 
wir Frau I. Engelfür die Ausführung einer großen Anzahl rechnerischer Kontrollen. 
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Berechnung des Wärmeüberganges an laminar umströmten 
Körpern mit konstanter und ortsveränderlicher 


Wandtemperatur* 
Von W. Dienemann in Stuttgart 
(Aus dem Institut für Strömungsmechanik der Technischen Hochschule Braunschweig.) 


Für die Berechnung des ebenen Problems der thermischen Grenzschicht an einem laminar umströmten 
Körper bei inkompressibler Strömung wird ein neues Näherungsverfahren angegeben, das in der rechne- 
rischen Durchführung wesentlich bequemer ist als die bisherigen Verfahren. Außer dem Fall des geheizten 
Körpers mit konstanter Wandtemperatur kann damit auch der Wärmeübergang an einem Körper mit orts- 
veränderlicher Wandtemperatur behandelt werden. Gleichzeitig ergibt sich auch ein vereinfachtes Be- 
rechnungsverfahren für die dreidimensionale Strömungsgrenzschicht an einem schräg angeblasenen Zylinder 
(schiebender Tragflügel). 

For the two-dimensional problem of the ihermal boundary layer in the case of laminar incompressible 
flow a new approximate method is given, which is much more convenient than the existing methods. Thus 
it is possible to treat the heat transfer on a body with constant or arbitrarity variable surface temperature. 
There results also a simplified method for the calculation of the three-dimensional boundary layer of the 
oblique flow round a cylinder (yawed wing). 

Pour le calcul du probleme plain de la couche limite thermique & un courant incompressible une nou- 
velle methode de rapprochement est indiquee, laquelle est beaucoup plus commode pour l’execution du cal- 
cul que la methode appliquee jusqu’a present. Outre le cas du corps chauffe & temperature constante de la 
paroi on peut aussi traiter la transition de chaleur a un corps & temperature non-stationnaire de la paroi & 
l’aide de cette methode. En mäme temps il s’ensuit aussi une methode simplifiee de calcul pour la couche 
limite iridimensionnelle du courani & un cylindre flanque de travers par le courant (aile portante eu 
derapage). 

Ina BEIINCHeHHA INIOCKOH 3amaym TePMuYecKOrO MOTPAHHYHOTO CIOA Y NAMHHAPHO oÖTe- 
KAEMOIO Tea HPH HEeCHKHMAECMOM TeYeHHH MAeTcH HOBEIH IIPHÖAMKeHHBIÄ CIIOCOÖ, KOTOPBLÜÄ, B 
OTHOIIeHHH pacYeTa, CYINECTBEHHO YAOÖHeeE IIPe>KHHX CITOCOÖOB. KpoMe cııyyaa HarpeBaeMoro 
Te1a C IOCTOAHHOÄ TeMNepartypoü CTeH, HO 3TOMY CIIOCOÖY MO:KeT TAKIKEe ÖBITb PACCMOTPeHa 
Tem1olepenaya Te1a C HECTAIHOHAPHOH TeMIepatypoä cTeH. Bmecte c TEeM NOAy4YaerTca 
TAK’Ke YIPOINeHHBIH CIIOCOO MIA BEIYHCHEHNA TPEeXMepHOrO IMOTPAHHYHOTO CIOA TeUeHHA y 
KOCO HATEKAEMO!O HMAHHAPA (ABHrAMIMTerO KPELIA). 


I. Einleitung 


Wird ein geheizter Körper von einer Flüssigkeit oder einem Gas umströmt, so gehen bei 
kleiner Zähigkeit und Wärmeleitfähigkeit des strömenden Mediums die Werte von Geschwindig- 
keit und Temperatur an der Wand in jeeiner dünnen Schicht in die Werte der Außenströmung 
über. Diese Übergangszonen mit großen Gradienten von Geschwindigkeit und Temperatur 
senkrecht zur Wand heißen Strömungs- bzw. Temperaturgrenzschicht. Die 
Berechnung des Wärmeüberganges von dem geheizten Körper in das strömende Medium er- 
fordert die Berechnung der Geschwindigkeits- und Temperaturverteilung in diesen Grenz- 
schichten. Ist die Strömungsgeschwindigkeit von der Größenordnung der Schallgeschwindig- 
keit, so daß das strömende Medium als kompressibel angesehen werden muß, so beeinflussen 
sich Strömungs- und Temperaturgrenzschicht gegenseitig. In diesem Fall wird die Wärme- 
bilanz außer durch Konvektion und Leitung maßgeblich durch die Reibungs- und Kompressions- 
wärme mitbestimmt. Überdies müssen in diesem Fall die Stoffbeiwerte Zähigkeit und Wärme- 
leitfähigkeit im allgemeinen noch als temperaturabhängig angesehen werden. Für kleine 
Strömungsgeschwindigkeiten dagegen (inkompressible Strömung, Mach-Zahl—0) und für 
mäßige Temperaturdifferenzen zwischen Wand und Außenströmung können Reibungs- und 
Kompressionswärme vernachlässigt und die Stoffbeiwerte als konstant angesehen werden. In 
diesem Fall wird die Strömungsgrenzschicht unabhängig von der Temperaturgrenzschicht, was 
die Berechnung erheblich vereinfacht. 

Ein Überblick über die Berechnungsverfahren der Strömungs- und Temperaturgrenz- 
schichten ist kürzlich von H. Schlichting [1] gegeben worden. Bei der Berechnung der 
Temperaturgrenzschicht haben sich die sog. exakten Lösungsverfahren, welche unmittelbar die 
Differentialgleichung für die Temperaturverteilung zu lösen suchen [2], als ebenso mühsam er- 
wiesen wie bei der Strömungsgrenzschicht. Man hat deswegen für die Temperaturgrenzschicht 
Näherungsverfahren von der gleichen Art ausgearbeitet, wie sie für die Strömungsgrenzschicht 
zuerst von K. Pohlhausen [3] angegeben wurden, und die sich dort als praktisch sehr 


*) Gekürzte Fassung der von der Fakultät Maschinenwesen der Technischen Hochschule Braunschweig 
enehmigten gleichnamigen Dissertation des Verfassers, 1951. Hauptberichter: Prof. Dr. H.Schlichting 
icberichter: Prof. Dr. E.Schmidt. Ein kurzer Auszug dieser Arbeit wurde bereits im Journal Aeronaut 
Sci. 18 (1951), S. 64 veröffentlicht. 
Y 
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brauchbar erwiesen haben. Ein solches Näherungsverfahren für den allgemeinsten Fall der 
kompressiblen Jaminaren Grenzschicht mit Wärmeübergang gab kürzlich J. Ginzel [4] an), 
während Näherungsverfahren für die Temperaturgrenzschicht bei inkompressibler Strömung 
schon früher von G. Kroujiline [5] und E. Eckert [6] mitgeteilt worden sind. Das 
Verfahren von Ginzel ist wegen der großen Allgemeinheit (kompressible Strömung, temperatur- 
abhängige Stoffbeiwerte) in der Durchführung ziemlich langwierig, aber auch die Verfahren von 
Kroujiline und Eckert für die Temperaturgrenzschicht in inkompressibler Strömung 
lassen sich noch wesentlich bequemer gestalten, wenn man die Vereinfachungen des Pohl- 
hause n-Verfahrens, de H.Holsteinund T. Bohlen [7] für die Strömungsgrenzschicht 
angegeben haben, in naheliegender Weise auch auf die Temperaturgrenzschicht anwendet. Die 
Ausarbeitung eines solchen gegenüber den bisherigen noch wesentlich vereinfachten Rechen- 
verfahrens für die Temperaturgrenzschicht bildet den wesentlichen Inhalt unserer nachstehenden 
Untersuchungen. 

Dabei wird sich auch gleichzeitig ein Beitrag zur Berechnung dreidimensionaler Strömungs- 
grenzschichten mitergeben. Bei inkompressibler Grenzschicht am schiebenden Zylinder ?) 
gilt nämlich für die Geschwindigkeitskomponente parallel zu der Erzeugenden die gleiche 
Differentialgleichung wie für die Temperaturverteilung des gleichen senkrecht angeström- 
ten Zylinders für die Prandtl-Zahl eins. Die von uns erhaltenen Vereinfachungen der 
Berechnung der Temperaturgrenzschicht ergeben deswegen auch gleichzeitig eine Vereinfachung 
der Grenzschichtrechnungen am schiebenden Zylinder, wie sie von W.R. Sears [8] und 
J.M. Wild [9] mitgeteilt wurden. 

Darüber hinaus konnte unser neues Näherungsverfahren noch auf einen Fall angewendet 
werden, der bisher erst wenig behandelt worden ist, aber doch großes praktisches Interesse besitzt, 
nämlich auf das Wärmeübergangsproblem eines Körpers mit örtlich veränderlicher Wandtempera- 
tur. In diesem Fall ergeben sich einige bemerkenswerte Besonderheiten für den Wärmeübergang: 
Für den lokalen Wärmestrom von der Wand in die Flüssigkeit oder umgekehrt gilt hier keineswegs 
mehr das Newtonsche Abkühlungsgesetz, welches diesen der Temperaturdifferenz zwischen 
Wand und Außenströmung proportional setzt. Es kann hier vielmehr an gewissen Stellen der 
Wand der Fall eintreten, daß dort die Wärme in der ‚‚falschen Richtung‘ strömt, d.h. daß die 
Wand z.B. Wärme aufnimmt, obgleich sie wärmer ist als die Außenströmung außerhalb der 
Reibungsschicht. Dies ist auf Besonderheiten der Temperaturverteilung in der Reibungsschicht 
zurückzuführen. Die Nusseltsche Wärmeübergangszahl, der ja das Newtonsche Ab- 
kühlungsgesetz zugrunde liegt, verliert hier völlig ihren Sinn. Ein erster Vorstoß zur Behand- 
lung dieser Vorgänge wurde von D.R. Chapman und M.W. Rubesin [10] gemacht, 
welche die längsangeströmte Platte in kompressibler Strömung auf Grund der Differential- 
gleichungen für die Temperaturverteilung exakt behandelten. Eine einfachere Lösung, ebenfalls 
für die Platte, wurde von H. 
Schlichting [11] angege- 
ben. Unser neues Näherungs- 
verfahren ist auf beliebige Kör- 
performen und beliebige Ver- 
teilung der Wandtemperatur 


Möglichkeit, dieses interessante 
Sondergebiet des Wärmeüber- 
ganges näher zu untersuchen. 


II. Die Grundgleichungen 


Wir betrachten die ebene 
laminare Strömung um einen 


Bild 1. Erläuterungsskizze. 


a = Grenze der Strömungsgrenzschicht, beliebig beheizten zylindrischen 
b = Grenze der Temperaturgrenzschicht, Ö 1 

ö = Dicke der Strömungsgrenzschicht, Körper, der senkrecht zu aeLueT 
ö; = Dicke der Temperaturgrenzschicht, Längsachse mit der Geschwin- 

U(x) = potentialtheoretische Geschwindigkeitsverteilung außerhalb der 1 : ö 1 
ee digkeit 1 angeströmt wird 
Tu@)= Temperaturverteilung an der Wand. (Bild 1). Es seien & und Y die 
rechtwinkligen Koordinaten 


längs der Körperoberfläche bzw. senkrecht dazu; der Staupunkt liege bi = y=0. Es sei 
T,,(x) die Wandtemperatur und 7, die konstante Außentemperatur in der ungestörten Strö- 


1) Ein analoges Näherungsverfahren für den Fall der wärmeundurchlässigen Wandist vonE.Grusch- 
witz veröffentlicht worden: Calcul approch6 de la couche limite laminaire en ö&coulement compressible sur 
une paroi non conductrice de la chaleur. ONERA publication No. 47 (1950). 

2) Zylinder, der nicht senkrecht, sondern schräg zu den Erzeugender angeströmt wird. 


anwendbar und gibt damit die, 


u a a; 
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mung; U(x) bedeute die Geschwindigkeit der Potentialströmung am äußeren Rand der Strö- 
mungsgrenzschicht. Es werden die folgenden Annahmen getroffen: 


1. Die Anströmungsgeschwindigkeit U, sei klein gegenüber der Schallgeschwindigkeit, so daß 
das strömende Medium als inkompressibel angesehen werden kann (Mach-Zahl Ma — 0). 

2. Die Zähigkeit u und die Wärmeleitzahl } seien klein als Voraussetzung für die Gültigkeit 
der Grenzschichttheorie. 

3. Der Temperaturunterschied zwischen Wand und Außenströmung T',(x) — T, sei überall 
klein gegenüber der absoluten Temperatur 7,, so daß die Stoffbeiwerte u und Aals kon- 
stant (temperaturunabhängig) angesehen werden können. Auch die Prandtl-Zahl 
Pr = v/a sei konstant. 

4. Reibungs- und Kompressionswärme werden vernachlässigt. 

Hiermit lauten die Strömungsgleichungen und die Energiegleichung mit den üblichen Grenz- 
schichtvereinfachungen, mit wund vals Geschwindigkeitskomponenten in der x- bzw. y-Richtung 
und mit T als Temperatur, (vgl. z. B. [1], Kap. 14): 


ou, MW er 
Sa (Kontinuität) EIER (1), 
2 
eo +v re + Da (Impuls «-Richtung) . . (2), 


u— +V— =a (Energiegleichung) . . . (3). 


Dabei ist a =A/o g c, die Temperaturleitfähigkeit, o die Dichte, g die Schwerebeschleunigung, 
v»— u/o die kinematische Zähigkeit und c, die spezifische Wärme bei konstantem Druck, bezogen 
auf die Gewichtseinheit. 


Die Randbedingungen sind: 
y-0:u=v=0; T=T,(@))| 
me = Un: T=J, Ya 


Die vom Körper in das strömende Medium übergehende Wärmemenge je Flächen- und Zeitein- 
heit ist nach dem Fourierschen Wärmeleitungstheorem 


oT 
ga(2) = HE RE N 
Dabei bedeutet (27/0y), den Temperaturgradienten in Richtung der Wandnormalen, genommen 
an der Wand. Es seig > 0, wenn Wärme vom Körper in die Flüssigkeit übergeht. 

Das Gleichungssystem (1) bis (4) ist so beschaffen, daß zunächst aus Gl. (1) und (2) die 
Strömungsgrenzschicht errechnet werden kann, die vom Temperaturfeld unabhängig ist. Nach- 
dem die Geschwindigkeitsverteilung in der Grenzschicht ermittelt worden ist, liefert Gl. (3) das 
Temperaturfeld und schließlich Gl. (5) den Wärmeübergang. 


(4). 


III. Näherungsverfahren für die Strömungsgrenzschicht 


Für die Berechnung der Strömungsgrenzschicht um einen beliebigen Körper (ebenes Problem) sind 
zahlreiche Näherungsverfahren angegeben worden, von denen sich dasjenige vonK. Pohlhausen [3] als 
besonders zweckmäßig erwiesen hat, vor allem nachdem es von H. Holsteinund T. Bohlen [7] in der 
rechnerischen Durchführung noch erheblich vereinfacht wordenist. Da dieses Näherungsverfahren 
im folgenden sinngemäß auch auf die Temperaturgrenzschicht übertragen werden soll, möge es hier kurz er- 
läutert werden. Der Grundgedanke dieses (und der meisten anderen) Näherungsverfahrens besteht darin, 
daß anstatt mit den Strömungsgleichungen (1) und (2) selbst mit einem daraus über die Grenzschichtbreite 
gewonnenen Mittelwert gerechnet wird. Dieser wird erhalten durch Integration von Gl. (2) über y von y = 0 
bis y= ®. Dies liefert in bekannter Weise den Impulssatz der Grenzschicht in der Form, (vgl.z. B.[1], Kap. 8): 

dd sc he im ; 
er CA Er Se er (6). 


Dabei bedeutet 92) die Impulsverlustdicke, ö*(x) die Verdrängungsdicke und 7,(z) die Wandschubspannung. 
Diese Größen hängen mit der Geschwindigkeitsverteilung u(z, y) in der Grenzschicht zusammen durch die 
Beziehungen: 
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Das Näherungsverfahren arbeitet nun so, daß für die Geschwindigkeitsverteilung u(z, y) ein plausibler Ansatz 
gemacht wird. Damit lassen sich die Größen ®, ö* und r, nach Gl. (7), (8), (9) ermitteln, und Gl. (6) liefert 
dann eine gewöhnliche Differentialgleichung für die Impulsverlustdicke #(x) (oder einen anderen geeignet ge- 
wählten Grenzschichtparameter). Bezüglich der Einzelheiten der Rechnung sei auf [1], Kap. 12, verwiesen. 
Nach K.Pohlhausen wird für die Geschwindigkeitsverteilung ein Polynom 4. Grades für den dimensions- 
losen Wandabstand n = y/ö angesetzt, wo ö(x) die Grenzschichtdicke der Strömungsgrenzschicht bedeutet: 


U 
Tr untw ta tun; O<Snsl. 


Für n> list w=U. Die Koeffizienten a,,. .., a, werden festgelegt durch gewisse Randbedingungen, die man 
für die Geschwindigkeitsverteilung vorschreibt. Im vorliegenden Fall sind es die vier Bedingungen: 
9°u dU ou Du 
NEE um y 0: U; Be u. ae are (10). 
Die Koeffizienten a,,...., a, lassen sich dann durch die Dimensionslose 
_.6?dU 
B93 : PÜn a nk a a cn ee (11) 
ausdrücken, so daß sich das obige Polynom für w/U in der Form schreiben läßt: 
# 
Em t+ANmM ..... 00 unsre een. 2) 
mit 
1 
Fn)=2n 2? +7; EAU er a a Ei) a (13). 


Die Geschwindigkeitsprofile in der Grenzschicht bilden also eine einparametrige Schar mit dem Form- 
parameter A. Durch Einführen von Gl. (12)in (7), (8), (9) ergibt sichnachHolstein-Bohlen[7Jaus 
Gl. (6) schließlich für die Größe 92v = Z die folgende gewöhnliche Diffferentialgleichung 


dZ F(x) 
ae ve 
mit 
=, 40203400 
BERTEaT DE  R e EEENE FTE ENE, (15) 


Die Funktion Fist eine universelle Funktion des Parameters x, der in einfacher Weise mit A zusammenhängt. 
Diese Funktion ist von der Körperform unabhängig, sieistin [1], Gl. (12, 35), angegeben. Die Anfangswerte 
im Staupunkt (U =0) ergeben sich aus der Bedingung, daß dort die Anfangsneigung (dZ/dx), einen endlichen 
Wert haben muß. Es muß deshalb dort # = Osein. Dies ergibt die Anfangswerte: 


Ao= 7,052; %0= 0,0770. 
Ein bequemes Integrationsverfahren für Gl. (14) findet man ebenfalls in [1]. 


IV. Das Wärmeübergangsproblem bei örtlich konstanter Wandtemperatur 
a)Das Näherungsverfahren 


Wir machen für diesen Abschnitt die Annahme, daß die Wandtemperatur längs der um- 
strömten Körperoberfläche konstant ist, 7,, = const. Ferner werde eine dimensionslose Tem- 
peraturverteilung £(x, y) eingeführt, indem die Übertemperatur 7T(x, y) — T, mit der vor- 
gegebenen konstanten Temperaturdifferenz T,, — T» dimensionslos gemacht wird: 


BEneT, 
Be er 


Die Integration der Energiegleichung (3) von y = 0 bis y = oo führt auf die zur Impulsgleichung 


der Strömungsgrenzschicht, Gl. (6), analoge ‚‚Wärmestromgleichung‘ der Temperaturgrenz- 
schicht: 


) a eiserne 


2 [na al), N 


Für die Temperaturverteilung wird das folgende Polynom in n, = y/ö, (6, = Dicke der Tem- 
peraturgrenzschicht) angesetzt: 


t=1—2n+2n—n; 0O<nm<1:.:..0:..0..00.0.llB) 
Es erfüllt die folgenden fünf Randbedingungen: 
0% ot 0% 
y=0: =, Dh Yyadı: b=ß; En im . 2419); 
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Die Temperaturprofile in verschiedenen Abständen & vom Staupunkt sind also nach Gl. (18) 
zueinander affin, im Gegensatz zu den Geschwindigkeitsprofilen nach Gl. (12). 

Führt man Gl. (18) in die Wärmestromgleichung (17) ein, so kann diese integriert werden. 
Bei den Näherungsverfahren von G. Kroujiline [5jJund N. Frössling [2]erhält man so 
eine sehr umständliche Differentialgleichung für das Grenzschichtdickenverhältnis 


A(®) = d,(@)/6(e) , 


vgl. E. Eckert [6]. Statt dessen führen wir analog zur Impulsverlustdicke, Gl. (7), eine 
„‚ Wärmeverlustdicke “ ! 


De... a „co 


ein. Für diese Größe ergibt sich dann eine wesentlich bequemer zu lösende Differentialgleichung 
als für A bei den bisher bekannten Näherungsverfahren. 

Die Ausrechnung der Wärmeverlustdicke mit den Ansätzen (12) und (18) für die Ge- 
schwindigkeits- und Temperaturverteilung ergibt: 


ee (*)anı = f IR (n) Od, + A [ En) td = HA) . . . (21). 
0 0 0 


Mit Gl. (20) und (21) folgt aus der Wärmestromgleichung (17): 
d, d 2 
da (U d,) ne Pr Hd: ’ 
oder 
ER 12) EEE 
2 de vde Pr 
Entsprechend der Holsteinschen Verbesserung des Pohlhausenschen Verfahrens für 
die Strömungsgrenzschicht, Gl. (15), führen wir jetzt noch ein 


DL eg (22). 


y 


und 


TUE ME SEHEN. 
rd de vda\d) \& 


Dann ergibt sich für © (x) aus (22) die gewöhnliche Differentialgleichung 


\-4:2.M N (24). 


% 


RO N ee ne (25). 
% U 


Hiermit haben wir für die Berechnung der Temperaturgrenzschicht eine Differentialgleichung 
erhalten, die ebenso einfach ist wie die Gl. (14) für die Strömungsgrenzschicht. Die hier auf- 
tretende Funktion F,(A, A, Pr) ist universell, d.h. unabhängig von der Körperform, ebenso 
wie die Funktion F (x) für die Strömungsgrenzschicht nach Gl. (14). Da die Prandtl-Zahlin 
F, sehr einfach eingeht, kann man sich damit begnügen, die Funktion H, (A, 4A) ein für allemal 
im voraus zu berechnen und F, jeweils aus der Beziehung 


ee 4 2 2 4 
F=7,4—24 A Hi=;, 
zu ermitteln. Die Funktion H,(A, A) erhält man durch Ausführung der Integrationen in 
Gl. (21) zu: 

1 


en 3 1 1 1 u: ) 


H, Te > Kr . . . . . . . . . (26) 


3) H; muß eine Funktion von A und 4 sein, weil 


EEE WR 
a a A ist, 
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| 21 | 1902 


4>1: HAN pt warnt 2... 07h) 
al N 
+Al507 1802: 7 81021 3094 2 
Wir setzen: 
H,=H,(2) FAR (A. oo. we ne 


Die hiernach berechnete Funktion H,- A ist in Bild 2 dargestellt, und in Tabelle 1 sind die 
Funktionen H,(A) und H,(A) angegeben. 


2,2 


2,0 


004 0,08 012 0,16 02 0,24 028 02 035 04 044 048 052 HM, 4 08 


Bild 2. Kurventafelder universellen Funktion Ht -Ain Abhängigkeit von A und A, nach Gl. (27a, b). 


A = Formparameter des Geschwindigkeitsprofils nach GI. (11) 
4= ö,/8 = Grenzschichtdickenverhältnis. 


Anfangswerte: Im Staupunkt (U =0) muß auch für die Temperaturgrenzschicht 
die Anfangsneigung dO/dx einen endlichen Wert haben. Dies erfordert in gleicher Weise wie bei 
der Strömungsgrenzschicht, daß im Staupunkt F, =0 


Tabellel. Dieuniversellen ; ; : : 
Funktionen HA mon ist. Damit folgt aus Gl. (16) die Bedingung 


für das Wärmeüberganes- ev 2 
problem bei ni Wand- Pr > AHA A) 22: (28), 
aa Bis KR IR NE ET a wobei der Anfangswert A, = 7,052 bereits aus der Strö- 
A H, 100- H, mungsgrenzschicht bekannt ist. Aus Gl. (28) ergeben sich 
0 0 { 0 Tabelle 2. Die Anfangswerte im Stau- 
02 0.0265 0.179 punkt 4, Hr,*%inAbhängigkeitvonder 
0.4 0.0521 0'286 Prandtl-Zahl Pr für das Wärmeüber- 
0,6 0.0761 NY 349 gangsproblem bei konstanter Wandtemperatur 
0,7 0,0873 0,356 
0.8 0.0980 0.363 Pr Ao Hio “to 
0,9 0,1080 0,365 Ten 
1,0 0,1175 0,364 ) oo 0,3333 © 
1,2 0,1345 0,353 0,350 2,0 0,2022 1,154 
1,4 0,1492 0,337 0,705 1,5 0,1789 0,5077 
1,6 0,1618 | 0,319 1,000 1,304 0,1670 0,3346 
1,8 0,1727 0,301 1,981 1,0 0,1432 0,1446 
2,0 0,1822 0,284 5,162 0,7 0,1123 0,0436 
2,5 0,2008 0,247 24,580 0,4 0,0722 0,00587 
3,0 0,2146 0,217 00 f) {) 0 
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die Anfangswerte des Grenzschichtdickenverhältnisses A, im Staupunkt in Abhängigkeit von 
der Prandtl-Zahl. Sie sind zusammen mit den sich daraus ergebenden Anfangswerten von 
H,(4As, 4) = H,, und x,(A,, A,) = x,, in Tabelle 2 angegeben. 


Durchführung der Rechnung: Für die Durchführung der Rechnung erweist 
es sich als zweckmäßig, dimensionslose Größen einzuführen. Die Längen werden bezogen auf 
eine Bezugslänge 1 (z.B. die Profiltiefe, Bild 1), die Geschwindigkeiten auf die Anströmungs- 
geschwindigkeit U,. Werden die dimensionslosen Größen mit einem * bezeichnet, so hat man 
Der Ur ZU/U., O0* = OU,/l. 

Gl. (25) wird in dimensionsloser Schreibweise: 

dO* F,(A,4, Pr) 
dgr U* 


Die Lösung der Differentialgleichung (25a) kann graphisch oder numerisch erfolgen. Der 


Anfangswert von ©* ist 
A I-dU 
0 ol | TFT 8: 
U, da)o 
* 


Die Steigung us der Integralkurve ©*(x*) ergibt sich nach folgendem Schema: 


. (25a). 


da* 
— u en FE nn en EEE 
Ze NR nm. au A| O* | n Pr dU 5 Mn En a H, = H, | 2x | nal a 
B ur: Ta | U Bild? | = a1.(26) | is 
R | | | | a, 
bekannt aus zeilenweise zu berechnen 


Strömungsgrenzschicht 


Ist nun ©*(x*) ermittelt worden, dann lassen sich die übrigen Grenzschichtgrößen leicht 
Y9* 
Ha 


De 
berechnen. Insbesondere ist die dimensionslose Wärmegrenzschichtdicke 7 YRe = 


wobei Re = = i die Reynoldssche Zahl bedeutet. 


Wärmeübergang: Die von dem umströmten Körper in die Flüssigkeit (oder 
umgekehrt) übergehende Wärmemenge pro Flächen- und Zeiteinheit beträgt an der Stelle «: 


oT 
1) = 1) =), 
Wird die örtliche Wärmeübergangszahl « dimensionslos gemacht, so erhält man die örtliche 
Nusseltsche Zahl nach Gl. (16) und (18): 
IE), ( ot ) 
Nu(z) = & een 


oder 


no Be er... (29). 

er ee, 

In bekannter Weise ist also die Nusselt sche Zahl umgekehrt proportional der Temperatur- 
grenzschichtdicke. 


b) Beispiele 

Das von uns im vorigen Abschnitt angegebene Näherungsverfahren soll jetzt mit exakten Lösungen, 
anderen Näherungslösungen und Messungen verglichen werden. 

1. Längsangeströmte ebenePlatte: Für dielängsangeströmte ebene Platte ist die exakte 
Lösung nach E.Pohlhausen [12] bekannt, vgl. auch [1], Kap. 14. Wir wollen den Vergleich der exakten 
Lösung mit der Näherung auf die Wärmeübergangszahlen beschränken. 

Nach der exaktenLösungist die örtliche Nusseltsche Zahl, bezogen auf den Abstand » von 
der Plattenvorderkante bis zur betrachteten Stelle x: 


NEE Aa) 1 ZB EEE Er Er (30) 


96 Dienemann, Berechnung des Wärmeüberganges an umströmten Körpern m Er Page 


wo Rey, — U„x/v, und der Beiwert a,(Pr) vonE.Pohlhausen[l2]angegeben worden ist, vgl. [1], S. 248. 


Es können die exakten Werte von a,(Pr) recht gut angenähert werden durch a, = 0,332 \Pr. 
Für de Näherungslösung erhält man die Dicke der. Strömungsgrenzschicht aus Gl. (14) zu 


37 2 
a1 Var van ie SEN 


—— Messung £.Schmidt 
u. K. Wenner [13] 
0,37 


-—-- exakte Rechnuny 
Frössling [2J 


0,367 
Nöherung E. Eckert 
[6] 
0,35 —:.—  Nöherung 
W Dienemann 
0,34 
Dar —— 
a Re = 170000 
0,32 x = 707000 
° = 64450 
o = 3980 
031 & = 15550 C 
230 
Re = Ba 
% 


Bild!3. Die örtliche Wärmeübergangszahl Nu (x) der längsange- 
strömten jebenen Platte mit konstanter Wandtem- 
peratrurin Abhängigkeit von der Prandtl-Zahl Pr. Ver- 
gleich der Näherung:lösungen mit der exakten Lösung nach Bild4. Örtliche Wärmeübergangszahl eines geheizten Kreiszylin- 
E. Pohlhausen [12], Näherung I nach Gl. (35); Näherung II ders; bei konstanter Wandtemperatur; Vergleich von 
nach GI. (53c). Theorie und Messung. Zylinderdurchmesser Z = 100 mm. 


Nöherung Eckert 

Näherung Kroujiline - Frössling 
Näherung Dienemann Abschn. 4 
Näherung Dienemann Abschn. 5 


identisch mit Kroyjiline - Frössling 


\ Eckert bis x*» 0,8 © 


Bild 5. Örtliche Wärmeübergangszahl Nu, Dicke der Temperaturgrenzschicht Ei und Grenz- 


schichtdickenverhältnis 4 an einem geheizten elliptischen Zylinder vom Achsenverhältnis 
a,/b,=4 beikonstanter Wandtemperatur. Bezugslänge für Re und Nu: 
U = große Achse deselliptischen Zylinders. 


Die Wärmeverlustdicke folgt aus Gl. (25) wegen 4 = O und A= 0 zu: 


oder 


2 02 
u. (*) ee: 


ee Fo 2 Zu De Zu SE amd » _ 
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Mit den Werten von H; nach G1.(27a, b) erhält man nach Division von G1.(32) durch GI. (31) die Bestimmungs- 
gleichungen für 4: 
1 


A < 1: 4%(1,135 — 0,182 4? + 0,047 48) = > . (33a), 
y 
1 1 1 
A . Auer ES 2 — a I ne . 
> 1:2,552 4? — 2,552 A + 1,135 — 0,1825 + 0,047, — 5- (33b) 
Die Lösung dieser beiden Gleichungen sei angesetzt in der Form - 
Alla OENB ee 
Für die Nusseltsche Zahl erhält man aus GI. (29) 
Nu(2) 2 0,343 
VRx 6YU,pe-4 Aus 
Bild 6. Geschwindigkeits- und Temperaturverteilung in der Grenz- 
(2) schicht eines geheizten elliptischen Zylinders vom Achsenyerhält- 
5 nis a/b, =4beikonstanterWandtemperatur fürdie 
U Prandtl-Zahlen Pr = 0,7; 1; 10. A = Ablösungspunkt. 


Us Lage der Profile an den folgenden Stellen: 

12 

z | Kr (8) = 
(1) (2) | (3) (4) | (5) (6) (7) | Ablösungs- 

10 | pkt. 

08 Pr A 0,285 | 0,536 0,899 


Bild 6a. Potentialtheoretische Geschwindigkeitsverteilung. 
In Bild 6e ist die Dicke ö der Strömungsgrenzschicht 


er — fach überhöht. 
100 ss] 


Ü 
es 1 
Fass 


0,6 


0,8 
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wenn man ö = 5,83» x2/U, nach GI. (31) einsetzt. Somit wird wegen Gl. (34): 
Nu(z) _ 0,343 
[Re  BiPr) 


In Bild 3ist die NusseltscheZahl über der Prandtl-Zahl aufgetragen, und zwar für die exakte Lösung 
nach G1.(30) und für die Näherung I nach Gl.(35). Die Übereinstimmung ist recht gut. 


2.Kreiszylinder: Als weiteres Beispiel für unser Verfahren wurde der quer angeströmte Kreis- 
zylinder gerechnet, weil in diesem Fall einerseits genaue Messungen ausgeführt worden sind [13], und weil 
andererseits hierfür auch die exakte Lösung von N. Frössling [2] und die Näherungslösung von 
E. Eckert [6] vorliegen. Ebenso wie Eckert haben wir bei unserer Rechnung die experimentelle Ge- 
schwindigkeitsverteilung zugrunde gelegt und die Prandtl-Zahl Pr = 0,7 (Luft) gewählt. In Bild 4 ist die 

aa) 1 
Ti er er 2 
Messungen verglichen mit unserem Näherungsverfahren, mit demjenigen von Eckert und mit der exakten 
Lösung vonFrössling. Die Übereinstimmung ist sehr gut. 

3. Elliptischer Zylinder: Schließlich haben wir unser Näherungsverfahren noch mit den 
Verfahren von E. Eekert [6] und Kroujiline-Frössling [2], [5] verglichen für den elliptischen 
Zylinder vom Achsenverhältnis 4: 1, der bei Pr = 0,7 parallel zur großen Achse angeströmt wird. Dabei ist 
die potentialtheoretische Druckverteilung und die Strömungsgrenzschicht nach H. Schliehting und 


VPr (Näherung DD... on 0. » 2 vr lab): 


örtliche Wärmeübergangszahl Nu = a: 1A = über dem Kreiswinkel aufgetragen. Dabei sind die 


BE 
d/l = 0,% 


’ 


07 0,2 0,3 0,4 0,5 0 02 0,4 0,6 0, 


|x 


Bild 7. Wärmeübergang in einem Schaufelgitter von symmetrischem Joukowsky-Profilen J 015 bei verschiedenem Teilungs- 
verhältnis £/2 und bei konstanter Wandtemperatur. 


a) örtliche Wärmeübergangszahl längs der Schaufeltiefe, Nu und Re bezogen auf !, Pr = 0,7. 
b) potentialtheoretische Geschwindigkeitsverteilung. 


A. Ulrich [14] zugrunde gelegt worden. Bild 5 zeigt den Verlauf der örtlichen Nusseltschen Zahl und 


der Verdrängungsdicke ö7 der Wärmegrenzschicht sowie auch das Grenzschichtdickenverhältnis A = ö,/ö). 
Bei der Wärmeübergangszahl und bei A ist die Übereinstimmung unseres Verfahrens mit den übrigen voll- 


kommen, bei öj ergeben sich für das Ec kertssche Verfahren gegenüber den anderen geringe Abweichungen, 
die aber ohne Bedeutung sind. Ein anschauliches Bild von der Geschwindigkeits- und Temperaturverteilung 
für dieses Beispiel gibt Bild 6, wo an verschiedenen Stellen längs der Kontur das Geschwindigkeitsprofil und 
die Temperaturprofile für die Prandtl-Zahlen Pr = 0,7; 1; 10 dargestellt sind. Damit die Profilformen 
unmittelbar verglichen werden können, ist neben der dimensionslosen Geschwindigkeitsverteilung u/U die 
dimensionslose Temperaturverteilung in der Form 1 —t = (T, — T)/(Tu — T„) dargestellt. Bei der längs- 
angeströmten Platte würden in dieser Darstellung für Pr =1 Geschwindigkeits- ne Temperaturprofil iden- 
tisch sein. Im vorliegenden Fall trifft dies nur näherungsweise zu, was auf den Einfluß des Druckgradienten 
auf die Ausbildung der Grenzschicht zurückzuführen ist. 


4. Schaufelgitter: Als weiteres Beispiel ist in Bild 7 der Wärmeübergang an den Schaufeln 
eines Gitters dargestellt. Das Gitter bssteht aus Schaufeln des symmetrischen Joukowsky-Profiles der Dicke 


&) Die Verdrängungsdicke der Wärmegrenzschicht ist definiert durch 
et 


ST To)dy= BEST 
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d/l= 0,15, es hat keine Staffelung und wird senkrecht zur Gitterfront angeströmt, so daß keine Ablenkung 
der Strömung vorhanden ist. Außer der alleinstehenden Schaufel (t/I— ») sind drei verschiedene Teilungen 
gerechnet worden. Obgleich die potentialtheoretische Geschwindigkeitsverteilung) sich mit der Teilung erheb- 
an ändert, ist der Einfluß der Teilung auf den Verlauf der örtlichen Nusseltschen Wärmeübergangszahl 
sehr gering. 

Auf Grund der Ergebnisse dieser Beispiele kann abschließend festgestellt werden, daß unser neues 
Näherungsverfahren zur Berechnung des Wärmeüberganges an laminar umströmten Körpern von beliebiger 
Form gute Ergebnisse liefert. Hinsichtlich der Genauigkeit ist unser Näherungsverfahren mit dem von 
Kroujiline-Frössling gleichwertig,in der rechnerischen Durchführung aber wesentlich einfacher. 
Beim Vergleich mit dem Eckertschen Verfahren ist das unsrige bezüglich der Genauigkeit ebenfalls als 
gleichwertig anzusehen. Da aber die von Eckert seinem Verfahren zugrunde gelegten Voraussetzungen nur 
als Hypothesen anzusprechen sind, die kaum zwingend begründet werden können, müssen die guten Ergeb- 
nisse seines Verfahrens als mehr oder weniger zufällig angesehen werden. Dies zeigt auch die Tatsache, daß 
bei seinem Verfahren bei der Strömungsgrenzschicht die Impulsverlustdicke und bei der Temperaturgrenz- 
schicht die Verdrängungsdicke die besseren Ergebnisse liefert, wie Eckert durch Probieren festgestellt 
hat. Demgegenüber beruht unser neues Verfahren auf dem physikalisch einleuchtenden Grundgedanken, 
daß die Gleichungen für den Impulsverlust der Strömungsgrenzschicht und den ‚Wärmeverlust‘‘ der 
Temperaturgrenzschicht erfüllt sein müssen, 

Ein ganz besonderer Vorteil unseres neuen Verfahrens ist jedoch, daß es auf allgemeinere Probleme der 
Temperaturgrenzschichten übertragen werden kann, die mit den bisherigen Näherungsverfahren nicht behandelt 
werden können. Dies ist insbesondere der Fall beim Wärmeübergangsproblem mit örtlich veränderlicher 
Wandtemperatur, das im nächsten Abschnitt behandelt werden soll. 


V. Das Wärmeübergangsproblem bei örtlich veränderlicher Wandtemperatur 
Das im vorigen Abschnitt für konstante Wandtemperatur angegebene Verfahren soll jetzt 
auf den Fall der örtlich veränderlichen Wandtemperatur erweitert werden. Die Wandtemperatur 
des Körpers sei eine beliebige stetige Funktion der Lauflänge x, also 7, (x). Die Temperatur- 
differenz gegenüber der Außenströmung werde angesetzt in der Form: 


BEER AT), After art (86), 


Dabei bedeutet AT,,o = Tao — Ts #0 die Temperaturdifferenz an einer beliebigen festen Stelle 
des Körpers; f (x)ist die dimensionslose Wandtemperatur. Der Gradient df/d« = f(x) sei überall 
vorhanden; er darf jedoch keine extrem großen Werte annehmen, damit die bei der Aufstellung 
der Energiegleichung (3) benutzte Annahme 9°?7/0x? < 0?T/dy? Gültigkeit behält. Dies ist er- 
füllt, wenn f/(z) — f (x)/ö, oder kleiner ist. 


a)Das Näherungsverfahren 
Wir führen eine dimensionslose Temperaturdifferenz gegenüber der Außenströmung analog 
zu Gl. (16) ein durch 


BET UT. To De 


Aamaıt, 7. TI. NT I u 


Für # lautet die der Gl. (17) entsprechende 
Wärmestromgleichung: 


d or’ 
|“ ey «(&), as) 
In Abschnitt IV wurden für die Tempera- 

turverteilung längs der Wand affine Profile ge- 
wählt, Gl. (18). Dies kann für konstante Wand- 
temperatur als recht gute Näherung angesehen 
werden, wenngleich es auch dort nach der 
exakten Lösung von N. Frössling nicht 
genau erfüllt ist. Bei örtlich veränderlicher 
Wandtemperatur dagegen ist die Affinität der 
Temperaturprofile auch in grober Näherung 
nicht mehr erfüllt. Dies läßt sich leicht ein- Bild 8. Temperaturgrenzschicht an einer längsangeströmten ebenen 
sehen am Beispiel der längs angeströmten late mit Sul ch vorn ihn hie 
ebenen Platte, deren Wandtemperatur vorn 2 Geschwindigkeitsverteilung in der Grenzschicht 
höher und weiter hinten tiefer ist als die 
der Außenströmung, Bild 8, vgl. auch [11]. Ganz vorn sind die Temperaturprofile vom Typus A, 
wie er bei konstanter Wandtemperatur geläufig ist. Die stromabwärts abnehmende Wand- 
temperatur, dT,„/dz< 0, beeinflußt zunächst nur die wandnahen Teile der Temperaturgrenz- 


5) Diese verdanke ich Herrn Dipl.-Ing. L. Speidel. 
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schicht. Infolgedessen ergibt sich an der Stelle C, wo die Wandtemperatur gerade mit der Außen- 
temperatur übereinstimmt, ein Temperaturprofil vom Typus C, dessen Gradient senkrecht zur 
Wand positiv ist, (07/dy), > 0, also das umgekehrte Vorzeichen wie an der Stelle A hat. Es 
fließt also hier Wärme von der Flüssigkeit in die Wand, obgleich die Wandtemperatur gleich der 
Außentemperatur ist. Daher muß es zwischen den Stellen A und C eine Stelle B geben, für die 
(oT/ay), = 0 ist, und für die somit der Wärmestrom nach Gl. (5) gleich Null wird, obgleich 
T, > Ts ist. Stromabwärts von der Stelle C ergeben sich Temperaturprofile mit negativer 
Temperatur in Wandnähe und positiver Temperatur weiter außen. Diese Überlegungen führen 
durch Vergleich mit den Verhältnissen in der Strömungsgrenzschicht mit Druckgefälle zu der 
Feststellung: 

Der Einfluß der örtlich veränderlichen Wandtemperatur auf die Temperaturverteilung in der 

thermischen Grenzschicht entspricht dem des Druckgradienten auf die Geschwindigkeitsvertei- 

lung in der Strömungsgrenzschicht. 

Deshalb muß für die Temperaturverteilung beim Wärmeübergangsproblem mit örtlich ver- 
änderlicher Wandtemperatur ein Ansatz mit einer einparametrigen Schar von Profilen gewählt 
werden wie beim Pohlhausen-Verfahren für die Strömungsgrenzschicht mit Druck- 
gradient, Abschnitt III. Dem Formparameter der Geschwindigkeitsprofile A nach Gl. (11) 
wird ein analog gebildeter Formparameter der Temperaturprofile entsprechen, welcher den 
Gradienten der Wandtemperatur längs der Wand, also f’(x), enthält. Wir setzen für die Tem- 
peraturverteilung ein Polynom 5. Grades an: 


Diihs 2 3 4 5 
DEE ’ = fa) t=f)otantamtemtamtan),; O<m<I.. . (89). 


Das Polynom +6: '''+c;n; hateinmal die gleichenfünf Randbedingungen zu erfüllen wie 
bei konstanter Wandtemperatur, Gl. (19). Darüber hinaus schreiben wir als sechste Bedingung 
noch die folgende Wandbindung vor: 


2a _ (®) AT, 
ee 0) 


die sich aus Gl. (3) durch Differentiation nach % ergibt, wenn man danach y = 0 setzt. Die Be- 
stimmung der Koeffizienten &, .. ., c, ergibt für die Temperaturverteilung 


12 =/)Fn)+AHm) - - - - (4) 
mit 
0,8 
Fı(n =(8 —5n+5m —3 m); | 
. (42). 
0,6 1 3 4 1} 
Gm.) TRRET: (m —6n +8 —3 nm) 
0,4 Dabei bedeutet 
5 A iD u(2+2)r@ . (43) ®) 
Pi BE 5 aM 
den Formparameter der Temperaturprofile. Die Funk- 
tionen F,(n,) und G,(n.) sind in Bild 9 dargestellt. 
Weiter wird nun wieder, wie in Gl. (20), die 
Wärmeverlustdicke 
Bild 9. ee * und @, für die j u 
rn = [(d-r)ay . . .. 2 004) 


- ö 
eingeführt, sowie auch H, = d,/ö, wie in Gl. (21). Damit ergibt sich aus der Wärmestrom- 


gleichung (38), wenn man noch a = Ei Ay] | ö, einführt: 


un -Ero+2]. 


6) Die Analogie von A; mit dem Pohlhausen-Parameter A nach Gl. (11) erkennt man noch 
deutlicher, wenn man 4; schreibt in der GEN 


nn a de (24 4a. 


a u a 


=. 
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Wird ferner 


DE (45) 
y 
gesetzt und 
MdU 
Ber m Gl u HAAICHIT N une a a (46 
leer 0.1 ia As H; (46), 
so erhält man für © (x) die gewöhnliche Differentialgleichung 
H, 10 A 
Ba 
ar Alla, AMAAPA|,.......«n. 


U U 


Hiermit haben wirfür die Berechnung der Temperaturgrenzschicht bei veränderlicher Wand- 
temperatur eine Differentialgleichung erhalten, die von der gleichen Art ist wie diejenige für 
konstante Wandtemperatur, Gl. (25). Die hier auftretende Funktion F,[f(®), A, A, A, Pr] ist 
universell, d.h. unabhängig von der Körperform, ebenso wie die entsprechende Funktion F, nach 
Gl. (25) bei konstanter Wandtemperatur. Da die Parameter f(x), Pr und A, in F, verhältnis- 
mäßig einfach eingehen, kann man sich damit begnügen, nur die Abhängigkeit der Funktion F, 
von den Parametern A und Aein für allemal im voraus zu berechnen, und F, jeweils aus der 
Beziehung 


_ Hrio Ir. 
ee am" 1 ame N len u |) A? hg A RI RTN 48 
nr —2 44H: (48) 


zu ermitteln. Dabei hängt H, außer von Aund A auch noch von f(x) und A, ab. 


Die Funktion H,[f(«), A, A, A] erhält man aus 


1 
DE Er 
ee, - v) an 
H, =f(«) [H,(A) + AB,(A)) + A [Hs (A) + AH (A) - - : » » - (49) 
LA. ir. (Aa). 


Die Ausrechnung mit «/U nach Gl. (12), (13) und /’ nach Gl. (41), (42) ergibt: 


10 ul u 1 


fee; 135 4° 
H,(d)= zelar 4-44 2-4 
Ardye zal1o54 f. TOR en ıeR A) 
a are l10s + - 4) ’ 


1 Del e IE I 
421: H,(d)= (4-m+37-nnt+ tea 


a iel DIR Uhl) 
ulesry (B-57 5047147 380 3) 
U au 9 1 2 1 Larl 
ul Aue 17, (16 "EI WATT 106%) 
are | or Mel mL 


ANA - — oz (7 - at aa) 


Die Funktionen H,,..., H, sind in Tabelle 3 angegeben. Die zusammengesetzten Funktionen Hi (4, d)= 
H,(4) + AH,(A) und H,(d, A)= H,(A) -+ AH,(A) sind in Bild 10 und 11 dargestellt. 


Anfangswerte: Im Staupunkt (U =0) muß wieder die Anlangsneigung do/da 
einen endlichen Wert haben. Dies erfordert, daß im Staupunkt F,= 0 ist. Der Formparameler 
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Tabelle3. DieuniversellenFunktionenH,H, der Temperaturprofile A, nach 


H, H, für das Wärmeübergangsproblem bei - ; “ 
” örtlich veränderlicher Wandtemperatur nach Gl. (50). Gl. (43) _ verschwindet im Stau 
punkt: A, =0. Der Formpara- 


A H, 12-H, | —ı®@-.A, |-10-H, meter der Geschwindigkeitsprofile 
hat den Wert A,= 7,052. Damit 

- 0 6 0 ß folgt aus Gl. (48) 

0,2 | 0,0315 0,210 | 0,042 0,264 1081 

0,4 | 0,0619 0,332 0,081 0,387 H,: f(x (>> —2 AH 4) 0 

0,6 | 0.0902 032 | 0118 0.418 J/a\FR - 

0,8 0,1158 0,411 0,149 | 0,396 

1,0 0,1384 0,406 | 0,174 0,349 oder 

1,2 0,1578 0,389 | 0,194 0,298 s 

1.4 0,1744 0,3677 | 0210 | 0,20 0,2363/[45 : Hıı (As, Ao)] = Pr 

1,6 0,1886 0,344 0,222 | 0,210 (51 

1,8 0,2006 0,322 0232 | - 0,178 2 ; 

2,0 0,2109 | 0,302 0,239 _ | 0,151 Hieraus ergeben sich die Anfangs- 

2 | an | re an 0.180 werte des Grenzschichtdickenver- 

6 | 0'253 0'253 | 0.098 hältnisses 4, (Pr) im Staupunkt. 

2,8 0.2402 0239 0,256 | 0.086 Sie sind zusammen mit den sich 

3,0 0,2456 | DO 0,259 | 0,076 daraus ergebenden Anfangswerten 


von H,(Av 40) = Ho und 
#.(Ao, 40) = %:0 in Tabelle 4 angegeben. Bemerkenswert ist, daß diese Anfangswerte von der 
Verteilung der Wandtemperatur f(x) unabhängig sind. © 


AHT| 


02 0,4 06 0,8 170 12 14 1,6 10 VZO 2 2% 


Bild 10. Die universelle Funktion H,,(4, 4) für das Wärmeübergangsproblem mit 
örtlich veränderlicherWandtemperatur nachGl. (49), (49a), (50). 


Durchführung der Rechnung: Es werden wie im Abschnitt IV dimensions- € 
lose Größen eingeführt: x* = a/l; U*= U/U,; 9*= OU,/l. Gl. (47) wird dann in dimensions- 
loser Schreibweise: 


d9* _ FıLf(ar), Au, A, 4, Pr] 


Tabelle4 Die Anfangswerteim Stau- da* U* 
punkt A,H:r, #0, inAbhängigkeitvonder 


Prandtl-Zahl Pr für das Wärmeüber. Die Lösung von (47a) kann wieder era- 
gangsproblem bei ortsveränderlicher Wand- 8 (22) 8 


temperatur phisch oder numerisch erfolgen. Der An- 
2 fangswert von O*ist 


. (47a). 


Pr As | His *to De 
| = old) 

0,254 2,0 0,2322 1,521 o bX/0 

0,327 1,8 0,2233 1,139 

a ne none a Die Steigung dO*/dx* der Integralkurve 

0.886 12 0.1852 0'348 O*(x*) muß etwas anders ermittelt wer- 

1.418 1.0 0,1666 0,196 den als bei konstanter Wandtemperatur, | 
2,552 0,8 0,1447 0,0945 da es jetzt nicht ohne weiteres möglich | 
vun nr ae ee ist, für irgendeinen beliebig gewählten 
39.383 0.3 0.0667 0.00282 Punkt x%*, 0* die zugehörige Steigung 
124.369 02 


0,0463 0,00062 dO*/dx* zu berechnen. Man muß vielmehr 


re aa Mach: Dienemann, Berechnung des Wärmeüberganges an umströmten Körpern 103 


bei einem vorgegebenen &* das Grenzschichtdickenverhältnis A annehmen und daraus den zu- 
gehörigen Wert von ©* ermitteln nach folgendem Schema: 


| 
se De y 
U Tr, A TV Re|f(=*) |f’(«*) 
— 
Rum mn m nn nn —{— — 
bekannt aus Strömungs- bekannt 
grenzschicht aus Vor- 
gegebener 
Wand- 
temperatur 
Be öt ——\® 
TYR=| A=(TYRe) 4 Hı = DR KR PA ELIA 
| I I- Hı,(4A, A)Hı,(A, A) nee Ne 4 * 
EN, A Bao | mia A re 
AT VRe x PrU* (? +) unse AH, ' U, das nach Gl. (47a) 


Vs 
zeilenweise berechnen 


Wärmeübergang: Wie D.R. Chapman und M.W. Rubesin [10] und 
H. Schlichting [11] ausführlich dargelegt haben, verliert die örtliche Nusseltsche 
Wärmeübergangszahl beiörtlich veränderlicher Wandtemperatur ihren Sinn, da die Wärmestrom- 


Az 10° 


wei 
— 
= 


02 


02 04 


6 08 10 12 74 7,6 18 20 32 24 


Bild 11. Die universelle Funktion H,,(4, 4A) für das Wärmeübergangsproblem mit 
örtlich veränderlicher Wandtemperatur nach GIl.(49), (49a), (50). 


dichte g(z) nicht mehr proportional der Temperaturdifferenz T,, — T,„ angesetzt werden kann. 
Man hat vielmehr g(z) aus Gl. (5) zu ermitteln und erhält im vorliegenden Fall mit Gl. (41) 
und (42) 


nn. .%% ..(82) 


A: ATyoYRe "yRe 


Ausdieser Gleichung liest man sofort ab, daß bei veränderlicher Wandtemperatur die Stelle, 
wo die Wärmestromdichte gleich Null ist, q = 0, nicht mit der Stelle zusammenfällt, wo die 
Wandtemperatur gleich der Außentemperatur ist, f(x) = 0. Vielmehr liegt die Stelle g(z) = 0 
bei in Strömungsrichtung abnehmender Wandtemperatur, f(x) < 0, A,<0, bei positivem f(x), 
d.h. dort, wo die Wand wärmer ist als die Außenströmung. Bei in Strömungsrichtung zunehmen- 
der Wandtemperatur [f/(x) > 0, A, > O]liegt die Stelle q(x) = 0 bei negativem f(x), also an einer 
- Stelle, wo die Wand kälter ist als die Außenströmung. — Diese Regel ist unabhängig vom Druck- 
gradienten der Außenströmung. _ 
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b) Beispiele 
1. Konstante Wandtemperatur 


Bei der Erprobung des neuen Näherungsverfahrens für ortsveränderliche Wandtemperatur mögen zu- 
nächst einige Fälle mit konstanter Wandtemperatur nachgerechnet werden, da man naturgemäß in diesem 
Fall auch befriedigende Ergebnisse erwarten muß. Bei konstanter Wandtemperatur ist wegen f’(z) = 0 nach 
G1.(43) auch Ar = 0. Dadurch vereinfacht sich das Formelsystem erheblich. Jedoch geht es nicht in das- 
jenige von Abschnitt IV über, da wir die Wandbindung Gl. (40), die sich jetzt auf (0°7'/dy?) = 0 reduziert, bei 
ortsveränderlicher Wandtemperatur neu eingeführt hatten. 

Die Formel für die örtliche Nusseltsche Zahl lautet jetzt in Abänderung von Gl. (29) ?) 


Nu _ 5/3 (52%) 
a e a). 
Y.Re —VRe 
3 Um Y ZH 
Buy Mar 


-15 


Bild 12. Temperaturfeld und Wärmeübergang an einer längsangeströmten ebenen Platte mit 
örtlich veränderlicher Wandtemperatur nach Gl. (54) für diePrandtl- 
Zahlen Pr = 0,7; 10. q(x) = örtlicher Wärmestrom, 


?) In Gl. (29) ist dabei die Zahl 2 durch 5/3 zu ersetzen, da die Wandtangente des TER DEE EEE 
nach (Gl. 41) und (42) jetzt 


oT\ __5Tu-To ., 
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Für das Beispiel derlängsangeströmtenebenenPlatte ergibt sich durch eine ähnliche 
Rechnung wie in Abschnitt IV, b1 als Bestimmungsgleichung für das Grenzschichtdickenverhältnis: 


A< 1: 4°(1,6226 — 0,2832 A? + 0,0756 49) = + N 2 cen.o., (53a) 
1 ZI 
A >1: 3,4083 4° — 3,0674 A + 1,1359 — 0,0945 40,83 up, - + : (53h), 
deren Lösung für verschiedene Prandt]- Zahlen wie in Gl. (34) in der Form 
A(Pr)= B(Pr) «Pr 
geschrieben wird. Damit ergibt sich die örtliche N usseltsche Wärmeübergangszahl aus Gl. (29) zu 
Nu(z) _ 0,2856 ,— 
ns) es ra. (680). 


YRez B(Pr) 


Das Ergebnis ist in Bild 3 als „Näherung II‘ miteingetragen; es stimmt mit der exakten Lösung etwas 
besser überein als die „Näherung I“ nach Abschnitt IV. 

Auch derelliptischeZylinder, der schon in Abschnitt IV gerechnet wurde, ist mit konstanter 
Wandtemperatur nach dem neuen Verfahren nochmals nachgerechnet worden. Die nach Gl. (52a) errechnete 
örtliche Nusseltsche Zahl ist in Bild 5 miteingetragen. Sie stimmt im ganzen Verlauf sehr gut mit der 
Lösung von Kroujiline-Frössling und mit der exakten Lösung überein. 


2. Ortsveränderliche Wandtemperatur 


Wir behandeln eine längsangeströmte Platte der Länge / mit einer Verteilung der Wand- 
temperatur nach der Gleichung: 


© 
Ty(z) ug Ir = ATyo » cos [" .) a u ee (54). 
Es ist hierbei Tyo— Ts = A To > 0 die Temperaturdifferenz zwischen Wand und Außenströmung an der 
Plattenvorderkante. Die Platte ist auf der vorderen Hälfte wärmer und auf der hinteren kälter als die Außen- 
strömung. Durch Vergleich mit Gl.(36) ist f(x) = cos (n + z/l). Die Rechnung wurde für die beidenPrandtl- 
Zahlen Pr = 0,7 und 10 ausgeführt. Das Ergebnis ist in Bild 12 dargestellt. Der qualitative Verlauf der 
Temperaturprofile wurde bereits auf 8. 99/100 besprochen. Die Übertemperatur der Wand gegenüber der Außen- 
Strömung geht bei x/! = 0,5 durch Null, dagegen der örtliche Wärmestrom an der Stelle z/! = 0,38. Zwischen 
diesen beiden Stellen nimmt die Wand Wärme aus der Flüssigkeit auf, obgleich sie hier wärmer ist als die 
Außenströmung. Für den gesamten Wärmestrom von der Platte in die Flüssigkeit, der sich durch Integration 
der q(z)-Kurve von x = O bis @ = lergibt, erhält man das unten angegebene Ergebnis: 


Werte von Q/A AT yo VRe: 


Pr gesamte Wärmemenge gesamte Wärmemenge Überschuß 
Platte — Flüssigkeit Flüssigkeit — Platte | Platte — Flüssigkeit 


— 0,544 


10 + 0,138 


0,7 | + 0,269 | — 0,242 | -+ 0,027 


-+ 0,682 


Obwohl also die Platte eine durch- 
schnittliche Wandtemperatur hat, 
die gleich der Außentemperatur ist, 
geht mehr Wärme von der Platte in 
die Flüssigkeit über als umgekehrt. 
Das Näherungsverfahren bie- 
tet auch die Möglichkeit, die Tem- 
peraturgrenzschichtt über z=|1 
hinaus in ein Gebiet konstanter 
Wandtemperatur weiter zu verfol- 
gen. Falls sich an die Temperatur- 
verteilung nach Bild 12 für #>1 
eine solche mit 7, — To = — ATyo 
anschließt, d.h. die Platte dort 
also kälter ist als die Außentempera- 
tur, so erhält man für Pr = 0,7 den 
Verlauf von g(x) nach Bild 13. Bei 
x=1|, wo sich die ortsveränder- 
liche Wandtemperatur an die kon- 
stante Wandtemperatur anschließt, 


ist eine ausgeprägte Spitze im Ört- Bild 13. Örtlicher Wärmestrom an einer längsangeströmten ebenen Platte mit vorn 
lichen Wärmestrom vorhanden, Zum örtlich veränderlicher Wandtemperatur und anschließender konstanter Wand- 


EN 
Vergleich ist in Bild 13 auch der ee 
örtliche Wärmestrom für von der Nase an konstante Wandtemperatur eingetragen worden. Demnach ist die 
Wirkung der ortsveränderlichen Wandtemperatur in der zehnfachen Entfernung nahezu abgeklungen. 

8 


——— konst Wandtemperatur 


Zee örtlich veränderliche 
Wandtemperatur 
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Um schließlich noch den allgemeinen Fall eines beliebigen Körpers mit ortsveränderlicher Wand- 
temperatur zu behandeln, ist der elliptische Zylinder von Abschnitt IV mit ortsveränderlicher 
Wandtemperatur berechnet worden. Die angenommene Wandtemperatur ist 


T 
Ty(z) -— To = ATyoo: 608 #7) a 
1 


Dabei bedeutet AT, > 0 die Übertemperatur im vorderen Staupunkt und /!, die Lauflänge vom vorderen 
Staupunkt bis zum laminaren Ablösungspunkt ®). Die Ergebnisse sind in Bild’ 14 dargestellt. Die Stelle, wo 
die Wärmestromdichte durch Null geht, liegt etwa am gleichen Ort wie bei der Platte. Wenn auch der Ver- 
lauf von q(z) um einiges von demjenigen an der Platte abweicht, so ändert sich an der Wärmemengenbilanz 


grundsätzlich nichts. Für die gesamte übergehende Wärmemenge Q/AAT'„. / Reı ergibt sich folgende Bilanz: 
Zylinder — Flüssigkeit: +0,331 
Flüssigkeit — Zylinder: —0,237 
Zylinder — Flüssigkeit: Überschuß: +0,094. 


Wiederum geht mehr Wärme vom Körper in die Flüssigkeit über als umgekehrt. 


48 


2,0 NR Mälds 

h A, Ei. SEREED AN 

Bild 15. Erläuterungsskizze. Strömung am 

ee schiebenden Zylinder. 

j ANzereeR 

BR . .. . 
ANLLITEN reimmmeente 

04 . u. . 


schrägangeströmten Zylinder _ 


Bei einem unendlich langen Zylinder, 
der nach Bild 15 schräg zu der Erzeugenden 


-94 angeströmt wird, wie es z.B. beim schie- 
benden Tragflügel der Fall ist, lauten die 
48 Grenzschichtgleichungen, vgl. [1], Kap. X, 
Bild 14. Temperaturgrenzschichtdicke & und örtlicher Wärme- au av 
a ae 4 4 nee tKonLiieatzer 
ratur nach Gl. (55). 0% 0Y 
ou ou dU ou 4 
- Impuls &-Richtun STD 
2 
u ec eg = ve (Impuls z-Richtung) . . . (58) 


mit den Randbedingungen: 
y=0 :uv=V=v=0 


y='ou = U: wi 


Dabei bedeutet w die Geschwindigkeitskomponente in der Grenzschicht in Richtung der Er- 
zeugenden (z-Richtung). Ferner ist die Annahme gemacht worden, daß die Außenströmung 
mit den Komponenten U und W von z unabhängig ist, so daß auch alle Grenzschichtgrößen 
von z unabhängig sind. 


°) Es ist I, = 0,8991, wo I die große Achse der Ellipse bedeutet. Bei dieser Wahl der Verteilung der 
Wandtemperatur lassen sich die Ergebnisse für den elliptischen Zylinder am besten mit denen der ebenen 
Platte vergleichen. 
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Wie man aus dem Vergleich der Gl.(1) bis (3) mit den Gl. (56) bis (58) ersieht, ist die Ver- 
teilung der w-Geschwindigkeitskomponente des schräg angeströmten Zylinders identisch mit 
der Temperaturverteilung des 


senkrecht angeströmten Zylin- h (LH 2 
ders im Falle 26 | ae 
L 10 
a Alasiik Pr =1. 2,4 | 9 
Damit kann man das in Ab- x 
schnitt IV für die Berechnung * N 7 
der Temperaturverteilung bei 
konstanter Wandtemperatur = x 
angegebene Rechenverfahren | 
sofort nutzbar machen für die 18 4 
Berechnung der Strömungs- Kagliien 
grenzschicht am schiebenden 16 a. hi 
Zylinder. Dadurch ergeben 
sich sehr erhebliche Verein- 14 Y 3 
fachungen gegenüber dem Nä- /} 
herungsverfahren von J.M. Y / -2 
wild [9 | M/X | 
Im einzelnen gestaltet sich 10 - /} -1 
die Rechnung folgendermaßen: DB ? 
Für die Geschwindigkeitskom- 0,8 | Je 25 PB AL 4:0 
ponente w gilt entsprechend % < 
Gl. (18) der Ansatz: 06 + BBur Senzuu 
a so BEE 


o<n,<1 


Nee ja jer 


AN\ EHEEEE 
IN 


BAER 


wobei n, = 4/6, ist, und 6,() 
die Grenzschichtdicke für w 
bedeutet, die von derjenigen 
für «, die Öö ist, verschieden 
sein kann. Zur Berechnung von 
ö, benötigt man die Impuls- 
gleichung für die z-Richtung, 
die man durch Integration von 
Gl. (58) erhält, in der Form: + Br 


e ri 
a Wwund ow ee 
73 u«W—u)dy=v Zeile 
“ 0) y Bild 16. Kurventafel der Funktion F,(A, 4,) nach G1. (63) zur Berechnung 
(60). der Strömungsgrenzschicht am schiebenden Zylinder, 


| 


N el 
N 
N 

DEKAN 
ZH 
ala 
ZA 

ZA, 

Ta AN, 

is 7 ir 
/EH 


N 
a 
I°7 


LINEEEBODNZBERSERER 
ia — ia 
EINER 


-04 


u w 
—-— | —I1— — ee nun, (61 
de: lit ve ar 
und analog zu Gl. (23) und (24): 
ds _y, oe ya, 
v An Wr da 


wobei A, = 6,/ö und H, = ®,,/ö, bedeuten. Man erhält dann als Differentialgleichung für Y(2): 


AA) |... ... (62). 
de U Si U 
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Dabei ist 
F (4,4) 22 H,0 A 0 Hy N 


identisch mit der universellen Funktion F, nach Gl. (26), wenn man dort Pr =1 setzt. Die 
Funktion F,(A, 4,) ist in Bild 16 dargestellt. Die Anfangswerte im Staupunkt (U = 0) folgen 


aus F.=0 zu: 
; A, = 1,3047 #0 20,102, = 0340, 
z 0 


vgl. auch Tabelle 2. 

Die Funktion H,(A, 4A,) ist identisch mit H,(A, A) nach Gl]. (27a) und (27b), wenn man A 
durch A, ersetzt. 

Da das Näherungsverfahren für den schiebenden Zylinder ein Spezialfall des Verfahrens 
von Abschnitt IV ist, bedarf es keines Vergleiches mit exakten Lösungen mehr. Das in Bild 6 
angegebene Beispiel für die Temperaturverteilung an einem elliptischen Zylinder für Pr =1 
gibt sofort die Geschwindigkeitsverteilung am schiebenden elliptischen Zylinder, wenn man in 


Bild 6 a = durch w/W ersetzt ?). 
10T ya ies 


VII. Zusammenfassung 


Für die Berechnung der Temperaturverteilung in der thermischen Grenzschicht an einem 
laminar umströmten Körper bei inkompressibler Strömung mit Vernachlässigung der Reibungs- 
und Kompressionswärme wird ein neues Näherungsverfahren nach v. Kärmän-Pohl- 
hausen-Holstein angegeben. Die Anwendung dieses Näherungsverfahrens auf das 
Wärmeübergangsproblem bei konstanter Wandtemperatur liefert Ergebnisse, die mit bekannten 
exakten Lösungen und anderen Näherungslösungen in guter Übereinstimmung sind, so daß es 
diesen Verfahren bezüglich der Genauigkeit als gleichwertig zu erachten ist. Bezüglich der 
Einfachheit der rechnerischen Durchführung, seiner physikalischen Begründung und der um- 
fassenderen Anwendung ist es den bisherigen Näherungsverfahren überlegen. Nach dem neuen 
Näherungsverfahren läßt sich auch das Wärmeübergangsproblem an umströmten Körpern be- 
liebiger Form mit ortsveränderlicher Wandtemperatur berechnen. Bei ortsveränderlicher Wand- 
temperatur ist es nicht mehr möglich, den örtlichen Wärmestrom durch die übliche Nusselt- 
sche Wärmeübergangszahl auszudrücken, da es hier u. U. Stellen gibt, wo Wärme von der 
Flüssigkeit auf die Wand übergeht, obgleich die Wandtemperatur größer ist als die der Außen- 
strömung. Auch für die Strömungsgrenzschicht an einem schiebenden Zylinder ergibt sich aus 
dem mitgeteilten Rechenverfahren eine beträchtliche Vereinfachung, da sich die Berechnung 
der Strömungsgrenzschicht am schiebenden Zylinder auf die Berechnung der Temperatur- 
verteilung am gerade angeströmten, geheizten Zylinder mit konstanter Wandtemperatur zurück- 
führen läßt, wenn dort die Prandtl-Zahl Pr =1 ist. 


Eingegangen am 24. März 1952. 
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BERINFEMITTEILUNGEN 


Ergänzungen zu dem Aufsatz ‚Der Energie- 
transport in strömenden Medien‘. 


In einem vorangegangenen Artikel!) hat der Ver- 
fasser zur Darstellung des mechanischen Energie- 
transportes, der in einem strömenden Medium statt- 
findet, einen Energietransportvektor eingeführt. Da 
dagegen von mehreren Seiten kritische Einwände er- 
hoben wurden, soll im folgenden eine direkt von den 
Grundgleichungen ausgehende Ableitung gegeben 
werden, die die Einführung eines solchen Vektors 
streng begründet. 

Wir gehen von den Stokes-Navierschen Glei- 
chungen aus und schreiben sie in der Form: ?) 


el + r (Ei)- ux(7 x n|=-eru+rr (1). 


Hierin ist o die Massendichte, v der Geschwindig- 
keitsvektor und U das Potential der Massenkräfte. 
Der Spannungsaffinor p hat die Form: 


=—-Jp+tnW,vp+u,P)—-JAVv.v (2. 


Hierin ist J der Einheitsaffinor, p der Druck, n und 
4 die dynamischen Zähigkeiten der Schub- und Kom- 
pressionsreibung. 

Um zu Energie- bzw. Leistungsgrößen überzugehen, 
multiplizieren wir Gl. (1) skalar mit vd und integrieren 
zeitfrei über ein beliebiges, festes Volumen des Strö- 
mungsraumes. 


Das Glied ov -[b x (7 x v)]ist identisch Null. 
Es gelten ferner folgende Identitäten: 


5 od _ © (olv?\ _|ol? de 
Hal 2 % 


er (er) Br en 


ov-VU=YV.(voeU)— UV-.(ov) 
v7 p=V (0 -pP)—(V -D)..p. 


Das Zeichen - - bedeutet das doppelt-skalare Produkt. 
Das Volumsintegral (3) gebt damit über in: 


er rre 


+7 .(eU)— UV:(ev) 


DH. p dr 


Bei Beachtung der Kontinuitätsgleichung: 


Bild 1 


Damit erhält man: 


Ne ter (5) er wxorx | 
. (3) 
+ov-7U—ov We 0 


dt bedeutet hierin das Volumselement. 
En Der Energietransport in strömenden Medien, diese Zeitschrift 
32 (1952), 9. 286. 

2). bedeutet Skalarprodukt, x Vektorprodukt, ‚dyadisches Pro- 
dukt. 


wird daraus, wenn man die reinen Divergenz-Gieder 
nach dem Gaußschen Satz in Oberflächen-Integrale 


umformt und au —= (0 annimmt: 


0]; 


fHettseo)-nf 


_2 [ (el! 
- 3 (eb! Heo)ar+ [uv,n.. pa: (8). 
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Hierbei ist do der nach innen gerichtete Vektor des 
Flächenelementes der Hülle. Rechnet man (A, v)..p 
nach Gl. (2) aus, so ergibt sich schließlich: 


[Heifso)-e 


° ((elvl? 
er & + e0)de- [pw war 


+[u, v)..nW, 0 +6, Milde — [aw-wrar (6). 


Die rechtsstehenden Glieder von Gl. (5) haben nun 
alle eine einfache physikslische Bedeutung und zwar 
der Reihe nach: die sekundliche Zunahme der Bewe- 
gungs- und der potentiellen Energie, die durch die 
Kompression bedingte thermodynamische Umsetzung 
und schließlich die Dissipation der mechanischen 
Energie durch Schub- und Kompressionsreibung. Im 
Sinne des Satzes von der Energieerhaltung kann man 
nun zwanglos das links stehende Oberflächenintegral 
als die durch die Hülle pro Zeiteinheit zuströmende 
mechanische Energie deuten. 

Der Energietransportvektor läßt sich somit durch 
den Ausdruck 


= u(ehl Heu) ae 


definieren und durch ein Flächenelement do wird 
während dt die Energie 


- (+ eu) war-do—p-d0)- var 2.22(8) 


transportiert. Wegen der Symmetrie des Affinors p 
gilt nämlich die Beziehung: 

(p-do)-v = (po) do 
Durch Einführen von Gl. (2) in Gl. (7) erhält © die 
Form: 


s-(2+l"4D)eu-nem0+nM 
2 Re as 


in Übereinstimmung mit G1.(8) der eingangs zitierten 
Abhandlung. 

Die durch ein Flächenelement hindurchtretende 
mechanische Energie setzt sich also nach Gl. (8) zu- 
sammen aus der Bewegungs- und der potentiellen 
Energie und der mechanischen Arbeit, die beim Durch- 
tritt durch doin der Zeit di von der durchströmenden 
Flüssigkeit nach außen geleistet wird. Dieser Ansatz 
wurde in der ersten Arbeit des Verfahrens an die 
Spitze gestellt und daraus eine Definition des Energie- 
transportvektors abgeleitet. Man erkennt überdies, 
daß die Definition von S nach G]. (7) für jedes Medium 
mit symmetrischem Spannungsaffinor gültig sein 
muß, 


Zum Abschluß sollen noch die Energietransport- 
linien für die Poisseuillesche Strömung einer in- 
kompressibeln Flüssigkeit durch ein zylindrisches Rohr 
dargestellt werden. Dabei soll eine Energietransport- 
linie als Einhüllende der Energietransportvektoren 
definiert werden. Es wurden folgende Annahmen zu 
Grunde gelegt: Wasser von 20° © strömt durch ein 
Kapillarrohr von « = 10cm Radius. Das Druck- 
gefälle J* = 3 liefert eine Strömung, deren Rey- 
noldssche Zahl Re = 75 beträgt. 

Am linken Ende von Abb.1 wird ein statischer 
Druck von 0,1 ata angenommen. Auf einer Strecke 
von ! = 34 cm ist dieser Druck auf 0 ata abgesunken. 
In Abb. 1 wurde der Radius 10®fach überhöht gezeich- 
net. Die Energietransportlinien kommen von der 
energiereichen Achsenzone und münden an der Wand, 
Dje Linien wurden auch im Zug-Teil der Strömung 


Bild 2 


strichliert gezeichnet, obwohl dieser Bereich physik- 
alisch nicht realisierbar ist. Man erkennt bereits aus 
Gl. (7a), daß der Transportvektor für große negative 
Werte von p gegen die Strömung gerichtet ist. Im 
Zugbereich strömt demnach die Energie in verkehrter 
Richtung, aber auch hier fließt sie von der Achse zur 
Wand. = 

In Abb. 2 ist die Umgebung des Überganges vom 
Druck- zum Zugbereich genauer herausgezeichnet. 
Hier wurde, bei gleichbleibenderÜberhöhung in radialer 
Richtung, eine Streckung um das 10*-fache in axialer 
Richtung vorgenommen. Wo die Summe aus Druck- 
und Bewegungsenergie verschwindet, besitzen die 
Transportlinien senkrecht zur Strömungsgeschwindig- 
keit gerichtete Tangenten, während sie im Druck- 
bereich spitze, im Zugbereich stumpfe Winkel mit ihr 
einschließen, 


Wien. G. Heinrich. 


BUCHBESPRECHUNGEN 


H. Hermes und H. Scholz, Mathematische 
Logik. (Enzyklopädie der math. Wissenschaften, 
Bd.I,, Heft 1, Teil 1.). 2. völlig neubearbeitete 
Aufl. 828. Leipzig 1952. B. G. Teubner Verlags- 
gesellschaft. Preis geh. 8,20 DM. 

InHöldersschönem Buche,, Die mathematische 
Methode‘ findet die Logistik eine zurückhaltende und 
z. T. ablehnende Beurteilung. Heute bildet die math. 
Logik einen Schwerpunkt der Diskussionen. Die 
lebhafte Entwicklung der Mathematik während der 
letzten Jahrzehnte hat dieses Forschungsgebiet be- 
sonders gefördert. Nachdem die Logistik ihre ersten 
praktischen Anwendungen gefunden hat, wird auch 
der praktische Mathematiker das Bedürfnis der 
Orientierung empfinden. 

In den gegenwärtigen Stand der math. Logik auf 
nur 82 Seiten einen verständlichen und umfassenden 


Einblick zu geben, war zweifellos eine außerordent- 
lich schwierige Aufgabe. Den großen Umfang an zu 
sichtender Literatur hat die math. Logik mit anderen 
Disziplinen gemeinsam. Hier mußte überdies in die 
der math. Logik eigene Symbolik eingeführt werden, 
und diese hat ihre besondere Geschichte. Ihre große 
Uneinheitlichkeit war früher ein oft geäußerter Grund 
der Ablehnung. 

Die vorliegende Darstellung ist eine recht glück- 
liche Lösung aller bezeichneten Schwierigkeiten, 
Auch der Nichtspezialist kann sich über jede ihn 
interessierende einschlägige Frage ohne weiteres 
orientieren. Die verhältnismäßig wenigen kalkül- 
mäßigen Voraussetzungen für das Verständnis ge- 
hören zur Eigenart der Mathematik überhaupt. Eine 
bis ins einzelne durchgeführte klare Aufgliederung 
des Stoffes und die dazu verwendete Technik leisten 
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bei der Lektüre vortreffliche Dienste. Wünschens- fahren und ihrer theoretischen Grundlagen. Es zeigt 


wert wäre es gewesen, auch auf die verschiedenen 
praktischen Anwendungen einzugehen, welche die 
mathematische Logik nisher bereits gefunden hat, 


Dresden, Draeger. 


Dr. K.H. Weise (o. Prof. a. d. Univ. Kiel), Vom 
mathematischen Denken. (Veröffent- 
lichungen der Schleswig-Holsteinischen Universitäts- 
gesellschaft. Neue Folge, Nr. 2.) 308. Kiel 1952. 
Verlag Ferd. Hirt. Preis geh. 1,20 DM. 

Bringt in gedrängter Kürze und klarer Gedanken- 
führung eine Oharakterisierung des math. Denkens. 
„Nicht allein das logische Schließen als math. Methode, 
sondern auch das Aufspüren von Gesetzmäßigkeiten, 
die in Grundannabmen niedergelegt werden, und die 
Untersuchung der Struktur des daraus hergestellten 
Gedankengebäudes machen das math. Denken aus.“ 
Das wird am Beispiel des Gruppenbegriffs, am Raum- 
problem und am Problem des Unendlichen erläutert. 

Der Vf. lehnt es einleitend ab, die ‚„„Bedeutung der 
Mathematik an ihrer Auswirkung in den Anwendungs- 
gebieten zu zeigen“. Gerade darauf kommt er nach 
der Klarlegung des Wesens math. Denkens zurück, 
gibt davon eine knappe, treffende Charakteristik und 
zieht zum Schluß die weitere Folgerung: ‚Hier liegt 
die Aufgabe vor, die Klarheit und vereinfachende 
Kraft modernen math. Denkens in geeigneter Aus- 
wahl zu vermitteln; denn es handelt sich dabei ja 
nicht um irgendeine weltfremde Gedankenspielerei, 
sondern um das Erkennen und Gestalten derjenigen 
Anlagen, die in dem heutigen naturwissenschaftlichen 
Zeitalter die geistige Orientierung des Menschen 
wesentlich mitbestimmen.‘ 


Dresden, Draeger. 


0. Haupt (o. Prof. a. d. Univ. Erlangen), Ein- 
führungin die Algebra undihre An- 
wendungen. (Mathematik und ihre Anwendun- 
genin Physik und Technik, Reihe A, Band 5, Teil 1.) 
2. Aufl, XVIII-+ 3708. Leipzig 1952. Akade- 
mische Verlagsgesellschaft. Geest u. Portig K. G. 
Preis geb. 23,— DM. 

Der spezifische Vorzug des Buches besteht darin, 
daß es dem Studierenden das Begreifen der ab- 
strakten Begriffsbildungen der modernen Algebra 
erleichtert, so daß er selbst, um mit dem Vf. zu 
sprechen, die Überzeugung gewinnt: „Die Abgren- 
zung, was konkret und was abstrakt, ist gefühls- 
mäßig bedingt, übrigens auch eine Funktion der 
Zeit“. An Bekanntes wird angeknüpft, und die 
Notwendigkeit der Entwicklung abstrakter Verall- 
gemeinerungen wird einsichtig. 

An allen wichtigen Stellen sind diesem Prozeß 
ausdrücklich prinzipielle Bemerkungen hinzugefügt. 
Darin geht der Vf. sogar erfreulicherweise über den 
rein algebraischen Bereich hinaus, wenn er sich z. B. 
auch über die axiomatische Methode, den Begriff der 
Existenz oder die Vollständigkeit von Axiomen- 
systemen äußert. Solche Hinweise sind wertvoll vor 
allem deshalb, weil sie hier aus dem Studium eines 
bestimmten mathematischen Stoffes erwachsen. 

Auch sonst gsschieht alles nur Mögliche zur Orien- 
tierung des Lesers, zur Lenkung der Lektüre und 
— durch Aufgaben, deren Lösungen am Schluß zu- 
sammengestellt sind — zur Festigung des jeweils 
Erarbeiteten sowie zur ständigen Selbstkontrolle. 


Dresden, Draeger. 


Dr.H.Bückner, DiepraktischeBehand- 
lung von Integralgleichungen. (Er- 
gebnis der angewandten Mathematik, Heft 1). VI+ 
127 S. mit 1 Abb. Berlin/Göttingen/Heidelberg 1952. 
Springer-Verlag. Preis brosch. 18,60 DM. 

Das Ziel der vorliegenden Monographie ist eine 
möglichst einheitliche Darstellung aller bisher bei 
Integralgleichungen 2. Art benutzter Lösungsver- 


sich, daß die meisten numerischen Methoden ent- 
weder Iterationsverfahren sind oder sich als Lösungen 
einer Integralgleichung mit Ersatzkern auffassen 
lassen. Auch die so wichtigen Variationsprinzipien 
führen — wenigstens formal — nicht aus diesem 
Rahmen hinaus. Mit der Reduktion der Aufgabe 
auf diese Kernfragen gelingt es, überall Zusammen- 
hänge zwischen bisher oft nur nebenher benutzten 
Formeln aufzuzeigen; an verschiedenen Stellen fin- 
den sich auch neue Ergebnisse. 

Nachdem im Abschnitt I die bekannten Formeln 
und Sätze aus der Theorie der Fredholmschen 
Integralgleichung zusammengestellt sind, werden im 
Abschnitt II Variationsprinzipien und Einschlie ßungs- 
sätze behandelt. Dieser Teil verdient besonderes 
Interesse, da er mit bisher unveröffentlichten Ar- 
beiten von H. Wielandt über Einschließungs- 
polynome den neuesten Stand der Forschung wider- 
spiegelt. 

Im III. Abschnitt werden die Iterationsverfahren 
in enger Anlehnung an bekannte Untersuchungen 
von Wielandt dargestellt. Weitere $$ beziehen 
sich auf die Behandlung des Spezialfalles hermi- 
tescher Kerne und Verallgemeinerung der klassischen 
Iterationsverfahren. 

Der IV. Abschnitt ist dem Verfabren mit Ersatz- 
kernen gewidmet. Nach allgemeinen Abschätzungs- 
sätzen für die Eigenwerte und Eigenfunktionen bei 
Kernänderungen werden zunächst Störungsverfahren 
behandelt. (Leider fehlt auch hier der Hinweis, daß 
mit einer Eigenfunktionsstörung jeweils zwei Eigen- 
wertschätzungen bestimmt werden können.) An- 
schließend wird die Approximation mittels ent- 
arteter Kerne untersucht. Hier ordnen sich auch 
die verschiedenen Rit zschen Verfahren ein. 

In einem letzten Abschnitt werden noch einige 
Integralgleichungen mit speziellen Kernen betrach- 
tet. 

Mit diesem umfassenden Inhalt empfiehlt sich das 
kurz und verständlich geschriebene Buch als z. Zt. 
modernste und umfassendste Darstellung der prak- 
tischen Methoden zur Behandlung von Integral- 
gleichungen 2. Art. Die Betonung liegt allerdings 
mehr auf der theoretischen als auf der numerischen 
Seite. 


Dresden. N.d. Lehmann. 


Dr. F. v. Krbek (Prof. a. d. Univ. Greifswald), 
Eingefangenes Unendlich. Bekenntnis 
zur Geschichte der Mathematik. VII + 331 S. mit 
128 Abb. Leipzig 1952. Akademische Verlagsgesell- 
schaft Geest & Portig K. G. Preis geb. 22,— DM. 

Das vorliegende Buchist KonradKnoppzum 
siebenzigsten Geburtstag gewidmet. Es behandelt 
in oft mehr oder weniger feuilletonistischer Art in 
einer Auswahl, die sich auf das Grundsätzliche be- 
schränkt, die Geschiche der Mathematik als eine 
Geschichte der Eroberung neuer Denkformen. Im 
ersten Abschnitt ‚Streifzüge in die Vergangenheit‘ 
schildert Verf. die Psychologie einer Anzahl von 
Forschern und ihre wechselvollen Schicksale. An- 
schließend folgt dann im Abschnitt ‚Pythagoras 
würde staunen‘ ein Überblick über die Entwicklung 
der Ideen der Arithmetik, Algebra und Analysis, und 
im dritten „Der überholte Euklid‘‘ werden Teile aus 
der Geschichte der Geometrie behandelt. Vor allem 
wird dabei die Bedeutung der Mengenlehre und der 
Gruppentheorie für viele Zweige der Mathematik 
aufgezeigt. Der Verf. hat nicht nur an den Fach- 
mann als Leser gedacht. „Einzelheiten wie Auf- 
fassung sind allgemein verständlich und für die breite 
Öffentlichkeit bestimmt“. Aber das Verständnis 
mancher Abschnitte fordert doch ziemliche Übung 
im mathematischen Denken und auch allerlei mathe- 
matische Kenntnisse. Für den Nichtfachmann wird 
dieLektüre noch dadurch erschwert,daßeine Reihe von 
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Versehen und Druckfehlern stehen geblieben sind, die 
der Mathematiker leicht richtig stellt. Z. B. ist in 
den Formeln auf S. 106 ein Vorzeichenfehler, auf 
Seite 134 wird festgestellt, „daß ein Zylinder einen 
dreimal so großen Inhalt besitzt wie die eingeschrie- 
bene Kugel‘, auf S. 181 fehlt ein Glied in der Formel 
usw. Einige Abb. sind sehr schlecht oder gar falsch, 
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wie Abb. 35, 72, 73, 78, 91, 96. Auf S. 172 wird auf 
eine Figur Bezug genommen, die offenbar nicht vor- 
handen ist, usw. Den Wert des Werkes als Ganzes 
beeinträchtigt das nicht. Man möchte dem äußerst 
anregend geschriebenen Buch viele Leser wünschen. 


Dresden. Willers. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


Dr. W. Reiechardt (Prof. a. d. Techn. Hochschule 


Dresden, Grundlagen der Elektro- 
akustik. VIII-+ 4648. mit 311 Abb. und 
mehreren Tabellen. Leipzig 1952. Akademische 


Verlagsgesellschaft Geest & Portig K. G. Preis geb. 
32,— DM. 


Dr. H. Hasse (Prof. a. d. Univ. Hamburg), 
Mathematikals Wissenschaft, Kunst 
und Macht. 34S. Wiesbaden 1952. Verlag für 
angewandte Wissenschaften. Preis geb. 4,80 DM. 


H. Reichardt u. W. Tollmien, Die Verteilung 
der Durchflußmenge in einem ebe- 
nen Verzweigungssystem. (Mitteilungen 
ausdem Max-Planck-Institut für Strömungsforschung 
Nr.7.) 488. mit 18 Abb. und 7 Tabellen. Selbst- 
verlag des Institutes. Göttingen 1952. Preis 5, — DM. 


Pref.Dr. U. Graf und Dr. #.J. Henning, Sta- 
tistische Metkoden bei textilen 
Untersuchungen. X 278 S. mit 71 Abb. 
Berlin-Göttingen-Heidelberg 1952. Springer-Verlag. 
Preis geb. 33,— DM. 


H.Drescher, Untersuchungenan einem 
symmetrischen Tragflügelmitspalt- 
los angeschlossenem Ruder bei ra- 
schen Änderungen des Rudereus- 
schlages (ebene Strömung). (Mittailun- 
gen aus dem Max-Planck-Institut für Strömungs- 
forschung Nr. 6). 71 8. mit 44 Abb. Selbstverlag des 
Institutes. Göttingen 1952. Preis 9,— DM. 


Prof. Dr. Fritz Neiß, Einführung in die 
Zahlentheorie. VIII -+113S. Leipzig 1952. 
S. Hirzel-Verlag. Preis kart. 4,80 DM. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


NACHRICHTEN 


Achter Internationaler Kongreß für Theoretische 
und Angewandte Mechanik. 


Vom 20. bis 28. August 1952 fand in den Räu- 
men der Universität Istanbul der Achte Inter- 
nationale Kongreß für Angewandte Mechanik 
statt. Für die deutschen Gelehrten war diese 
Veranstaltung auch deshalb so bedeutungsvoll, 
weil es der erste internationale Mechanikkongre3 
seit 1938 war, an dem deutsche Wissenschaftler 
wieder teilnehmen konnten. Etwa 470 Gelehrte aus 
30 Ländern waren der Einladung der türkischen 
Fachgenossen gefolgt, und es wurde von allen Teil- 
nehmern sehr schmerzlich empfunden, daß Gelehrte 
aus einer Reihe von Ländern, deren Beteiligung 
nicht minder erwünscht gewesen wäre, dieser so über- 
aus harmonisch verlaufenen Zusammenkunft fern- 
geblieben waren. 

Die rund 240 Vorträge waren auf Sektionen, bis 
zu 7 an der Zahl, aufgeteilt, und das ungewöhnliche 
Faktum ist hervorzuheben, daß man von dem An- 
hören eines Vortrages in einer Sektion unmittelbar 
in eine andere Sektion überwechseln konnte, weil 
der Zeitplan genauestens mit Hilfe von Weckeruhren 
eingehalten wurde. Außer den Sektionsvorträgen 
fand an jedem der 7 Sitzungstage ein allgemeiner 
Vortrag (‚General Lecture‘‘) zu einer von anderen 
Vorträgen freien Zeit statt, dessen Besuch also allen 
Teilnehmern möglich war. Als Vertreter des deut- 
schen Kulturkreises hielt der Unterzeichnete einen 
allgemeinen Vortrag über Schwingungen in laminaren 
Strömungen und die Theorie der Turbulenz. Dankens- 
werterweise ist geplant, die Verhandlungen des Kon- 
gresses bald im Druck herauszubringen, wobei die 
„General Lectures‘ sowie die ‚„‚Sectional Addresses‘“ 
in vollem Wortlaut wiedergegeben werden sollen, 
während für die Sektionsvorträge Auszüge vorgesehen 
sind. 

Näch der vielseitigen und anstrengenden wissen- 
schaftlichen Tagesarbeit vereinten sich die Mitglieder 
des Kongresses zu den freien Zeiten auf Empfängen 


und sonstigen gesellschaftlichen Veranstaltungen. Es 
war erhebend, zu beobachten, wie nicht nur die tür- 
kischen Kollegen sich die größte Mühe um die Or- 
ganisation des Kongresses gegeben hatten, sondern 
daß die ganze ehrwürdige Weltstadt Istanbul sich als 
Gastgeberin des Kongresses fühlte, ja man möchte 
fast sagen, die gesamte Türkei; denn der Kultus- 
minister gab einen Empfang für die Kongreßteil- 
nehmer, und an einem anderen Abend führten Volks- 
gruppen aus allen Teilen des Landes ihre Tänze vor. 

Die deutschen Teilnehmer wurden von allen Seiten 
mit großer Herzlichkeit wieder offiziell in die inter- 
nationale Gemeinschaft der Mechanikforscher aufge- 
nommen. Der Unterzeichnete wurde einstimmig indas 
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und Angewandte Mechanik gewählt. Die deutsche 
Sprachewar neben Englisch, Französisch und Italie- 
nisch eine der offiziellen Kongreßsprachen und wurde 
nicht nur von Mitgliedern, deren Muttersprachesie ist, 
sondernauch von Vertretern andererVölker, wie Skan- 
dinaviern und Türken, in ihren Vorträgen verwandt. 

Auf die wissenschaftlichen Ergebnisse des Kon- 
gresses soll erst eingegangen werden, nachdem die 
Verhandlungsberichte gedruckt vorliegen. Hier 
möge nur noch einmal der herzliche Dank an die 
türkischen Gastgeber zum Ausdruck kommen. 

Im Anschluß an den Kongreß wurden von der 
Naturwissenschaftlichen Fakultät der Universität 
Istanbul die Professoren Richard v. Mises, 
John v. Neumann, Ludwig Prandtl 
und Sir Geoffrey Taylor in feierlicher Zere- 
monie zu Ehrendoktoren ernannt. 


Göttingen. W. Tollmien. 


Würzburg: DozentDr.HerbertBilharz, 
bisher Universität Freiburg i. Br, wurde vom 
Bayerischen Staatsminister für Unterricht und 
Kultus zum planmäßigen außerordentlichen Pro- 
fessor für Mathematik an der Universität Würzburg 
ernannt, 
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